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Abstract 

This review deals with the theory of four-fermion interactions in curved spacetime. Starting with 
the D-dimensional Minkowski spacetime (2 < D < 4) the effective potential in the leading order of 
1/A-expansion is calculated and the phase structure of the theory is investigated. Using the same 
technique the effective potential for composite operator tpt/j in four-fermion models is calculated 
under the following circumstances: 

a) D-dimensional weakly curved spacetime (in linear curvature approximation), 

b) D-dimensional de Sitter and anti-de Sitter universe, 

c) Z)-dimensional Einstein universe. 

The phase structure of the theory is investigated analytically as well as numerically. Curvature 
induced phase transitions are discussed where fermion masses are dynamically generated. 

As an extension of four-fermion models we consider the gauged Nambu-Jona-Lasinio (NJL) 
model, higher derivative NJL model and supersymmetric NJL model in weakly curved spacetime 
where the effective potential is analytically evaluated. The phase structure of the models is again 
analyzed and the condition for the chiral symmetry breaking in the gauged NJL model is given in 
an analytical form. 

Finally the influence of two external effects (non-zero temperature and gravitational field, non- 
trivial topology and gravitational field as well as magnetic and gravitational field) to the phase 
structure of four-fermion models is analyzed. The possibility of curvature and temperature-induced 
or curvature- and topology-induced phase transitions is discussed. It is also argued that the chiral 
symmetry broken by a weak magnetic field may be restored due to the presence of gravitational 
field. Some applications of four-fermion models in quantum gravity are also briefly investigated. 

This paper was published in Prog. Theor. Phys. Suppl. 127 (1997) 93. 
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10 Summary and outlook 



1 Introduction 



An idea of dynamical symmetry breaking (DSB) was introduced by Nambu and Jona-Lasinio in quantum 
field theory in Ref. [pj (for modern reviews see, e.g., Ref. pi). The original Nambu- Jona-Lasinio (NJL) 
model and its generalizations (all of which are somehow related to the theory of superconductivity a 
la Bardeen, Cooper and Schriffer are extremely useful toy models in the study of composite states 
in various circumstances. To discuss the dynamical symmetry breaking nonperturbative treatments are 
essential. One of the simplest nonperturbative techniques is the 1 /./V-expansion scheme which works well 
in the four-fermion models. 

In the simplest version of the NJL model jij (or four-fermion model) the chiral symmetry of the theory 
is spontaneously broken according to the emergence of a vacuum condensate of the composite field of 
fermions tp '■ i' l P' t P} 0- At the same time the non-vanishing fermion mass is generated. In the spectrum 
one finds a mode of the composite scalar field : a ~ tptp, which may correspond to the Higgs field in the 
theory with the elementary Higgs scalar field. Clearly the model may be regarded as a composite Higgs 
model based on the mechanism of dynamical symmetry breaking. 

Four-fermion models (gauged four-fermion models) have been often considered as an effective low 
energy theory for quantum chromo dynamics (QCD, see Ref. Q for an introduction). The chiral sym- 
metry of such a theory is broken by the creation of quark-anti-quark condensates and quark masses are 
dynamically generated. Such a theory is extremely useful in dealing with non-perturvative phenomena 
of QCD. 

Moreover the (gauged) four-fermion models may be well applied to the description of the dynamical 
symmetry breaking in standard model (SM). At the electroweak scale the dynamical symmetry breaking 
may be caused by the fundamental fermion condensate as in the technicolor model. || The similar 
mechanism could work in the symmetry breaking at the GUT era in early universe. 

The description of elementary particle theories (QCD, SM, GUT, • • •) in terms of the effective four- 
fermion model is particularly convenient when the theories are considered under some external conditions 
(finite temperature, finite density, non-trivial topology, external gauge and gravitational fields, etc.). 
As an extreme external condition one immediately considers the early universe where the gravity is 
strong and the temperature is high. Since the dynamical symmetry breaking mechanism replaces the 
spontaneous symmetry breaking caused by the elementary Higgs scalar field, it is natural for us to 
employ the composite field er as an inflaton in the inflationary universe. Through renormalization 
group arguments one finds the necessity of the scalar-gravitational non-minimal coupling £ [0] whose 
value is governed by the renormalization group running. The value of £ varies in different theories and 
is crucial to have a successful inflation. 

One realizes an importance of studying the phase structure of the four-fermion model in curved 
background and finds the cosmological motivation to deal with this model. 

The main purpose of the present article is to present what is and is not established in the study of 
the four-fermion models in curved spacetime in full detail. We do not consider in the present article 
any cosmological applications. Instead we give a review of the phase structure of the different kinds of 
the four-fermion model with different gravitational backgrounds. Hence the classical external gravity is 
considered here as a probe of quantum four-fermion models and is regarded as a source of the curvature- 
induced phase transition in these models. 

The paper is organized as following. In §2 fundamental properties of four-fermion models are briefly 
reviewed. To investigate the phase structure of the four-fermion theory the effective potential is intro- 
duced in curved spacetime. We adopt the 1/N expansion method as a nonperturbative approach and 
calculate the effective potential in the leading order of the 1/N expansion. As an instructive rehearsal we 
work in the Minkowski spacetime and evaluate the effective potential analytically as well as numerically. 
On the basis of the effective potential we discuss the phase structure of the four-fermion model in the 
Minkowski spacetime. 

In §3 the weakly curved spacetime is considered. We keep only terms independent of and linear in the 
curvature. Evaluating the effective potential in the leading order of 1/N expansion the phase structure 
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of the theory is clarified and the curvature-induced phase transition is discussed in arbitrary dimensions 
(2 < D < 4). 

In §4 we consider the spacetimes with the constant curvature (de Sitter background, anti-de Sitter 
background and Einstein universe background) . In these spacetimes we can evaluate the effective poten- 
tial of the four-fermion model without making any approximation with regard to the curvature. With 
the aim of justifying the validity of the weak curvature approximation we investigate the phase structure 
in the constant curvature spacetime. We then compare the result with the one in the weakly curved 
spacetimc. 

Section 5 is devoted to the study of the gauged NJL model, supersymmctric NJL model and higher 
derivative four-fermion model in the weakly curved spacetime. We evaluate the effective potential in the 
gauged NJL model with finite cut-off and let the cut-off tend to infinity. In order to do so we use the 
renormalization group method and appeal to the equivalence with the gauge Higgs- Yukawa model. The 
analytical form of the condition for the chiral symmetry breaking is found. For the four-fermion model 
with higher derivative terms we also study the effective potential and the phase structure of the model. 
The phase structure of supersymmctric NJL model is also discussed in detail. 

In §6 the combined effect of curvature and temperature is discussed in curved spacetime. To introduce 
temperature we suppose the thermal equilibrium and we restrict ourselves to the Einstein universe 
background. The effective potential is evaluated at finite temperature in spacetimes with positive or 
negative curvature. For sufficiently high temperature we show that the broken chiral symmetry is restored 
and the symmetric phase can be realized even in the spacetime with negative curvature. 

In §7 we discuss the influence of the non-trivial topology to the chiral symmetry breaking in the weakly 
curved spacetime. We start with the calculation of the effective potential in the four-fermion model in 
the flat torus-compactified spacetime. We find the same quantity as in the weakly curved spacetime with 
the torus-compactified dimension. The possibility of curvature and topology-induced phase transition is 
explicitly shown. 

Section 8 is devoted to the investigation of chiral symmetry breaking in the weakly curved spacetime 
under the influence of the external magnetic field. The effective potential is found and its phase structure 
is studied. The chiral symmetry is shown to be broken in the NJL model under the external magnetic 
field. However, the inclusion of the external gravitational field with positive curvature acts in opposite 
direction and it may restore chiral symmetry. 

In §9 some applications of four-fermion models in quantum gravity are presented. In particulary, we 
discuss semiclassical approach in 2D Gross-Neveu-dilaton model. The conformal factor dynamics based 
on four-fermion theory is studied in 4D quantum gravity. The effective potential for composite gravitino 
in N = 1 supergravity on de Sitter background is also calculated. 

In §10 we give an outlook and prospects for future researches. 

Our notation is basically in conformity with the (— , — , — ) convention in the book by Misncr, Thornc 
and Wheeler. § 

2 Four-fermion models in Minkowski spacetime 

The four-fermion model is one of the prototype models of the dynamical symmetry breaking. In the 
present section wc briefly review the fundamental properties of the four-fermion models. A four-fermion 
model and its effective potential are introduced in curved spacetime. We adopt the l/N expansion as 
a non-perturbative approach to investigate the dynamical symmetry breaking in the model. The phase 
structure of the theory is discussed in Minkowski spacetime as an instructive example of our method. 
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2.1 Gross-Neveu type model 

Here we consider the Gross-Neveu type model in curved spacetime and discuss the symmetry and the 
symmetry breaking in arbitrary dimensions (2 < D < 4). It is defined by the action 



-gd D x 



N - A ( N 

k=l \fc=l 



(1) 



where index k represents the flavors of the fermion field "0, N is the number of fermion species, g the 
determinant of the metric tensor g llUl 7^ the Dirac matrix in curved spacetime and VpV' the covariant 
derivative of the fermion field tp. In two spacetime dimensions the theory is nothing but the Gross-Neveu 
model, j^] For simplicity, we neglect the flavor index below. 
The action (Q) is invariant under the discrete transformation 

W> — ► , / 2 \ 

07m i> — " V^m^ ■ 

In two or four spacetime dimensions this transformation is realized by the discrete chiral transformation 

ip — > 750 • (3) 

Below we call this Z 2 symmetry the discrete chiral symmetry. 

The discrete chiral symmetry prohibits the fermion mass term. If the composite operator constructed 
by the fermion and anti-fcrmion develops the non-vanishing vacuum expectation value 

W)^0, (4) 

fermion mass term appears in the four-fermion interaction term and the chiral symmetry is broken down 
dynamically. 

We consider the theory with N flavors of the fermion fields. Thus the theory also has SU(N) flavor 
symmetry 

tp — ► e 4 ^-'f 9aT ip , (5) 

where T a are generators of the SU(N) symmetry. Under the circumstance of the global SU(N) flavor 
symmetry we may work in a scheme of the 1/N expansion. 

For practical calculations in four-fermion theories it is more convenient to introduce auxiliary field 
cr.|| The generating functional is given by 

Z = J [d^}[dij}e' s , (6) 

where we choose the path-integral measure in terms of the field variables 



(7) 

-gy, 

so that the covariance under the general coordinate transformation fl(i|| is guaranteed. 
We consider the Gaussian integral 



C = j [da'} expz J ^-gd D x f-^^) 



(8) 



Obviously C is a constant. The path integral measure is invariant under the redefinition of the parameter 
a' 

a> ^a = a' (9) 



G 



If 



. A 

1 N 



After the redefinition, Eq. (||) becomes 

C = J [da]expi j v -.</'/"■'■ 



(10) 



As Eq. ([To]) is a constant, we are free to insert it in the right-hand side of Eq. (Q) without any change. 
We obtain 

Z = \ /[#][#] [d*]e iS «, (11) 



where action S v is given by 



N 



2An 



(12) 



The four-fermion interaction term in Eq. (^) has been canceled out by introducing the Gaussian integral 
©. 

The action S y is equivalent to the original action (|l|).[pT[ If the non- vanishing vacuum expectation 
value is assigned to the auxiliary field a 

(a)=m^0, (13) 

there appears a mass term for the fermion field tp and the discrete chiral symmetry (the Zi symmetry) 
is eventually broken. 

ip propagator: 
a propagator: 
iptpa vertex: 



ip loop: 

In Fig. 1 we present the Feynman rules for the action ( |l2| ) in Minkowski spacetime. These Feynman 
rules are useful to make quantum calculations in the Gross-Neveu type model. The Feynman rule for 
the a propagator contains the factor 1/N. Thus the radiative corrections including the a propagator are 
suppressed at the large N limit. Internal a lines have no contributions in the leading order of the 1/N 



7 



Figure 1: Feynman rules of Gross-Neveu type model in Minkowski spacetime. 



expansion. Therefore we easily calculate the radiative correction in the Gross-Neveu type model at the 
leading order of the 1/N expansion. For example the effective potential in Gross-Neveu type model is 
calculated from the two-point function of massive free fermions (see §2.2). 



2.2 Effective potential 

We would like to find a ground state of the system described by the four-fermion model. For this 
purpose we evaluate an effective potential for the field a. The ground state is determined by observing 
the minimum of the effective potential. If the auxiliary field a develops the non-vanishing vacuum 
expectation value, the ground state is no longer invariant under the discrete chiral transformation and 
then the chiral symmetry is broken down dynamically. In the present subsection we briefly review how 
to calculate the effective potential in general situation, i.e., curved spacetime. The case of the Minkowski 
space is included as a special case. 

We start with the generating functional 



Z[J] = / [#][<#] [da] exp iS y + iN / y/=y~d D xa{x) J(x) , 



(14) 



where J(x) is the source function for the field a{x). Performing the integration over the fermion fields 
tp and ip, we get 



Z[J] = J [da] exp iN J y^gd D x + ^j) - iln Det(i 7 Al V M 

We divide the field a into two parts 



a = ao + a . 



(15) 



(16) 



where ao is a classical background which satisfies the classical equation of motion and a is a quantum 
fluctuation. In terms of ao and a Eq. (|l5|) is rewritten as 



Z[J] = exp iN 



J J~gd D x (-^- + <r jj - iln Det(*y - a ) 



I [da] exp iN 
—iln Dct 



2A 



a 



(i7"V M - a ) 



(17) 



As may be seen through Fig. 1 radiative corrections in a are higher order corrections in the 1/N 
expansion. In the leading order of the 1/N expansion the generating functional W[J] for connected 
Green functions is given by 



NW[J] 



-ilnZ[J] 



8 



-0 



-gd D x 
1 

N 



2A 



a J - tin Det(r)^V M - <r ) 



(18) 



We pull out an obvious factor iV in defining the generating functional W[J]. 
The effective action T[a c } is defined to be the Legendre transform of W[J] 

T[a c ] = W[J] - J y/=gd D x cr c (x)J(x) , 

where the new variable a c is given by 

SW[J] _ / 1 

<7 c = — -!— - = a + O — 
5 J \N 



(19) 



(20) 



which corresponds to the vacuum expectation value of a in the presence of the source J. We substitute 
Eqs. (ph and @ into Eq. (|1 to get 



r[j] 



-gd D x—^- - tin Det(i7^V M - cr c 
2Ao 



The effective potential V(ct) (with TV factored out) is defined by 

I» 



V(a) 



(21) 



(22) 



where we put a c (x) = a, a constant independent of the coordinate x. The effective potential gives the 
energy density under the constant background a. The vacuum state should minimize the effective po- 
tential. Thus the effective potential is useful to determine the vacuum state in a static and homogeneous 
universe where the assumption a c (x) = a, a constant, is adopted. According to the Schwinger proper 
time method 113] the second term of the right-hand side of Eq. (21) for constant a reads 



In Dct(i7 M V M - a) 

= tr J d D x\n(ij tl V t , - a) 

= — tr J d D x\J—g J dsS(x, x; s) + const , 



(23) 



where const is a constant number which is independent of the field a and S(x, x; s) is the spinor two-point 
function which satisfies the Dirac equation 



(ifV p -«)S(i,jfi«) 



-S°(x,y). 



(24) 



Thus the effective potential (|22|) reads 

1 



V(a) = 



2A n 



-a — itx 



dsS{x, x;s) + O 



N 



(25) 



It should be noted that the effective potential is normalized so that V(0) = 0. From Eq. ( p5j ) we 
recognize that the effective potential is described by the two-point function S(x, x; s) of the free fermion 
with mass s. 
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We turn our attention to the theory governed by the action 



f N \ 2 / N 



~2N 



\k=l 



\k=l 



(26) 



In four dimensions the action (26) describes the NJL (Nambu-Jona-Lasinio) model. jlj The action is 
invariant under the chiral U(l) transformation in even dimensions, 



fl>. 



(27) 



The chiral U(l) symmetry prevents the action from having mass terms. Using the auxiliary field method 
and applying the similar method mentioned above, we obtain the effective potential for the NJL model: 



itr 



dsS(x, x;s)+0[ — 



(28) 



(29) 



We need, corresponding to the two four-fermion interaction terms in the action (|2q) , two kinds of auxil- 
iary fields a and w. If a' develops the non- vanishing vacuum expectation value the chiral U(l) symmetry 
is broken down and then Nambu-Goldstone (NG) mode appears |ll| (i.e., tt becomes masslessQ). In 
two dimensions it is well-known that the it loop has an infrared divergence and we cannot neglect the 
next-to-leading order terms of the 1/N expansion. Jl4| Therefore our present analysis is not sufficient for 
dealing with the action j26[ ) in two dimensions. 

The effective potentials ( pq ) an d y2g ) have the same form in the leading order of the 1/N expansion. 
Thus, using the effective potential (25), we can discuss the phase structure of both models except in two 
dimensions. 



2.3 Effective potential in Minkowski spacetime 

An instructive example to study the phase structure of the four-fermion model may be found in an 
infinite volume flat spacetime (i.e., Minkowski spacetime). In the present subsection we calculate the 
effective potential V(a) and show the phase structure of the Gross-Neveu type model in Minkowski 
spacetime. 

The two-point function S(x, x\ s) in Minkowski spacetime is given by 

S(x,x;s) = J ^ D]i l _ s , (30) 



where \J—g = 1. Inserting Eq. (|30|) into Eq. Q25J) and performing the integration, we obtain the effective 
potential in Minkowski spacetime, 

where trl is the trace of the unit Dirac matrix, the suffix for Vo(tr) is introduced to keep the memory 
that R = and L — > oo with R the spacetime curvature and L the size of a space or time component. 



4 In Gross-Neveu type model the symmetry under the transformation (gj) is not a continuous symmetry. Thus no NG 
boson appears. Our analysis is valid for Gross-Neveu type model in two dimensions. 
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The effective potential ( |3l| ) is divergent in two and four dimensions. It happens to be finite in three 
dimensions in the leading order of the 1/N expansion. As is well-known, the four-fermion model is 



renormalizable in the two-dimensional Minkowski spacetime. Therefore the potential (31) is made finite 
at D = 2 by the usual renormalization procedure. For D = 3, four-fermion model is known to be 
renormalizable in the sense of the 1/iV expansion. Jill The effective potential ( pl| ) happens to be finite 
in three dimensions in the leading order of the 1/N expansion.^] In four dimensions the four-fermion 
model is not renormalizable and the finite effective potential cannot be defined. In the present section 
we regard the effective potential for D = 4 — e with e sufficiently small positive as a regularization of the 
one in four dimensions and consider the theory for D = 4 - e as a low energy effective theory stemming 
from more fundamental theories. 

We introduce the renormalization procedure by imposing the renormalization condition, 



d 2 V (a) 



da 2 



A 



(32) 



where fi is the renormalization scale. From the renormalization condition ( p2[) we get the renormalized 
coupling A, 



1 

Ao 



A 



(4tt 



trl . . f D\ , 



(33) 



Replacing the bare coupling constant Ao with the renormalized one A we obtain the renormalized effective 
potential 



V (a) 



1 

2A 



^r 2 ,,D-2 

u n 



trl 



2(4tt) d / 2 
trl / D 

r 1 

(4 7r )«/2£» ^ 2 



(D-1)T 1 



D 
~2 



(34) 



In Minkowski spacetime the renormalized effective potential Vo(cr) is no longer divergent in the whole 
range of the spacetime dimensions considered here, i.e. for 2 < D < 4. 

We consider the two-, three- and four-dimensional limit of the effective potential (|^). Taking the 
two dimensional limit, D —> 2, we get 



_ 1 (ZY trl 
/i 2 2A \n J 8tt 



Taking the three dimensional limit D —> 3, we find 



V D=3 



1 

2A 



trl 

4-7T 



-3 + In 



5 1 1 

3 



(35) 



(36) 



Equations J3q) and (Q) are not divergent for finite a confirming that the divergent part in the effective 
potential (pW is removed by the renormalization procedure for 2 < D < 4. 

If we take the four-dimensional limit D — > 4, the effective potential (M) reads 



D=4 



2A 



a 



div 



trl 

4(4tt) 



e a 



div 



2 Vm 




(37) 



5 If we perform the momentum integration_by using the cut-off regularization, linear divergence appears in three di- 
mensions. The divergence is dropped in Eq. ( |3l| ) by the analytic continuation of T function. After renormalization both 
regularization methods give the same results. 
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where we express the divergent part by 

^div = 



4- D 



- 7 + In 47r + 1 



(38) 



On the other hand we can calculate the effective potential for D — 4 by using the cut-off regularization. 
Inserting Eq. (pOl) into Eq. (p5j) at D = 4 we obtain the effective potential for D = 4, 



1 



itr 



d 4 /c 



1 



(2tt) 4 # 



We cut off the higher momentum region of the divergent integral in Eq. (B9[) and get 



Vr 



to 



2A 



trl 

4(4tt) 



n Jl2 x 4 
Z<7 A (J 

2 



cr 4 In 



A 2 y 



(39) 



(40) 



where A is a cut-off scale of the divergent integral. We apply the same renormalization condition^ as 
shown in Eq. (p2|) and obtain the renormalized coupling A: 



1 

Ao" 



EL 

A 



trl 



(41) 



Substituting Eq. (|l]) into Eq. 



j) we obtain the renormalized effective potential for D = 4 



V D=Ai 




(42) 



We find that there is a correspondence (l6) between Eq. ( |37| ) and Eq. (|||) if we make a replacement 



C, 



div 



In 



(43) 



This correspondence rule is to be modified in curved spacetime (see §3.2). 



2.4 Phase structure in Minkowski spacetime 

In terms of the effective potential Vo(a) we can argue the phase structure of the four-fermion model. 
Using Eq. ( |34|) the effective potential is calculated numerically. In Fig. 2 we present the typical behavior 
of the effective potential evaluated at D — 2.5. In drawing Fig. 2 we introduce, for convenience, the 
new variables mo which will be defined by Eq. ( |47| ) and m' Q defined by 



trir(l -D/2) A 



i-D+i 



l/0?-2) 



(44) 



As is shown in Fig. 2 the shape of the effective potential is of a single and a double well for sufficiently 
small and large coupling constant A respectively. Hence we expect that the effective potential changes 
its shape at a critical value of the coupling constant Acr- We can see in Fig. 2 that the ground state 

6 Since the four-fermion model is not renormalizable in four dimensions, we do not apply the renormalization in four 
dimensions. Note here that the direct comparison of the dimensional reguralization with the cut-off rcguralization is 
possible only after renormalizing the coupling constant A under the common renormalization condition. 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

cr/m' 

^ < Acr 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

a/m 

A > Acr 



Figure 2: Behavior of the effective potential in Minkowski spacetime for D = 2.5. 



is invariant under the discrete chiral transformation for a sufficiently small coupling constant. On the 
other hand for A > Acr the minimum of the effective potential is located at non-vanishing a. In the 
latter case the auxiliary field a develops the non-vanishing vacuum expectation value and the discrete 
chiral symmetry is broken down dynamically. As the result the fermion acquires the dynamical mass. 

Next we calculate the dynamically generated fermion mass for A > Acr- It is equal to the vacuum 
expectation value of a which is obtained by observing the minimum of the effective potential. The 
necessary condition for the minimum is given by the gap equation: 



dV (a) 



0(7 



m /(m , A) =0, 



(45) 



where /(mp, A) is given by 



/(m ,A) = -ii 



D-2 



trl 



a (D-i)r(i 



D 



.D-2 



trl 



(47r) B / 2 



D 



r 1- — m 



D-2 



(46) 



The dynamical fermion mass is obtained by the non-trivial solution of the gap equation (45). Differen- 
tiating the effective potential ( |34| ) and solving the equation /(mo, A) = for mo, we find the dynamical 
fermion mass: 



m = fi 



(Air) 1 



trir(l - D/2) A 



-, l/(D-2) 



(47) 



In Fig. 3 we plot the dynamical fermion mass m as a function of the coupling constant A for 
D = 2, 2.5, 3, 3.5. The dynamical fermion mass decreases as 1/A increases, so that the dynamical fermion 
mass disappears at a critical value of A. Since no mass gap is observed at the critical point in Fig. 3, 
the phase transition at the critical coupling Acr is of the second order. In two dimensions the dynamical 
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Figure 3: Dynamically generated fermion mass for D = 2, 2.5, 3, 3.5. 



fermion mass decreases exponentially as 1/A increases. Thus the critical coupling constant Acr = for 
D = 2. 

For second order phase transition the critical coupling constant Acr is obtained by taking the limit 
ma — ► 0: 

lim /(m , Acr) -0, (48) 

mo— *-0 

where /(mo, Acr) is defined in Eq. (pq). Solving Eq. ( [48| ) we easily find the critical coupling constant 

i -i 



(4vr) D / 
trl 



(I- LhT[l-^ 



(49) 



For A > Acr the dynamical fermion mass is generated and the discrete chiral symmetry is broken down. 
In Fig. 4 we plot the critical coupling constant Acr as a function of the spacetime dimension D. 
For some specific values of D the solutions too and Acr simplify: 



m = ^e 1 - 2 -/^ 1 *), A 



Taking the four-dimensional limit, D — > 4, we find 



cr = ; D = 2, 

Acr-^-; D = 3. 
trl 



(50) 



a 



div 



In 



mo 



1 

Acr 



toq 
fj, 

3trl 

(4^)2 



trlA 



C div q 



3 a 



div 



where the divergent part Cjjy is defined in Eq. (j38|). 



(51) 
(52) 
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Figure 4: Critical coupling constant as a function of the spacetime dimension D. 



2.5 Renormalization group analysis in Minkowski spacetime 

We analyze the /3-function for the four-fermion coupling A. It has been often pointed out that the 
fermions acquire the dynamical mass for a larger coupling than the ultraviolet stable fixed point. S, [lTj 
The renormalization group /3-function for the four-fermion coupling is defined by 



Taking into account Eq. ( |33| ) with Eq. ( [49|) we obtain 

D - 2 



/3(A) = 



A(Acr — A) . 



(53) 



(54) 



in the leading order of the 1 /N expansion. 

In Fig. 5 we show the behavior of the /3-functions. As is obviously observed in Fig. 5 the theory 
is asymptotically free for D = 2 and the theory for 2 < D < 4 has a nontrivial ultraviolet stable fixed 
point at 

A = Acr ■ 



It simplifies for some special values of D: 

/3(A) = 



trl 

~ 2n 



A 



Taking the four-dimensional limit, D — > 4, we find 

3trl 



D = 2, 
D = 3. 



(55) 



(56) 



/3(A) = -2A 



(4tt) 2 \ C div 



(57) 
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Figure 5: Behavior of the /^-function for D = 2, 2.5, 3, 3.5. 



where C^ v 1S defined in Eq. (p8[). For a larger coupling than the ultraviolet stable fixed point a tachyon 
pole appears in the a propagator. Then the classical vacuum state becomes unstable. The generation 
of the dynamical fermion mass helps recovering the stability of the vacuum. 

The basic properties of the four-fermion models have been briefly reviewed. We calculated the effective 
potential and discussed the phase structure of the four-fermion model (Gross-Neveu type model) in 
Minkowski spacetime. In the leading order of the 1 /N expansion the effective potential is given by the 
integration of the spinor two-point function S(x,x;s) over mass variable s. Evaluating the effective 
potential in the leading order of the 1/N expansion it was shown that the discrete chiral symmetry is 
broken down dynamically for a sufficiently large coupling constant A > Acr- Dynamically generated 
fermion mass mo and the critical point Acr was obtained analytically by solving the gap equation ( f45|) . 
Since no mass gap appears at the critical point, the phase transition from the broken phase to the 
symmetric phase is of the second order. Applying the renormalization group analysis it was shown that 
the critical coupling constant Acr corresponds to the ultraviolet stable fixed point. 

In the present section we basically used the dimensional reguralization method to investigate the 
phase structure for 2 < D < 4 and considered the theory for D = 4 — e with e sufficiently small positive 
as a regularized version of the one in four dimensions. As was shown the effective potential can be 
calculated also by cut-off regularization. In this case the same result is obtained with recourse to the 
relation between e and A ([43]). 

In the following sections we will apply the similar method to the case of the curved spacetime. 

3 Four-fermion models in weakly curved spacetime 

At an early stage of the universe it is generally assumed that broken symmetries are restored and 
the physical phenomena are described by a more fundamental theory with a higher symmetry. The 
spontaneous symmetry breaking may be induced under the influence of the strong curvature, finite 
volume, non-trivial topology and so on. 

Here we discuss the curvature induced phase transition in four-fermion models. For this purpose the 
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effective potential^] of the four-fermion models need to be calculated in a general curved spacetime. In 
the present section we consider the theory in weakly curved spacetime. We shall confine ourselves to the 
case that the spacetime is curved very moderately and shall neglect terms involving derivatives of the 
metric higher than second order. 

There are many studies of the phase transition in weakly curved spacetime. Spontaneous symmetry 
breaking of the massless scalar theory with At/) 4 interaction has been studied in four dimensions (see for 
example Refs. |l|, and the references there in). In the scalar theory the continuous U(l) symmetry 
is restored for a sufficiently large positive curvature. The pioneering works in the four-fermion model 
have been performed by Itoyama,[^o) and Buchbinder and Kirillova plf in two dimensions. In Refs. 
J2^1~[pi| four-fermion models have been studied in three, four and arbitrary dimensions (3 < D < 4) 
respectively. In these literatures it is found that the possibility of the curvature-induced phase transition 
in four-fermion models. 

In this section we mainly follow Refs. |^| and 0] and discuss the phase structure of the Gross-Neveu 
type model in weakly curved spacetime by using the similar method explained in the previous section. 



3.1 Riemann normal coordinate expansion of S(x,y;s) 

As is shown in the previous section the effective potential of the Gross-Neveu type model is described by 
the spinor two-point function S(x,x; s) in the leading order of the 1/iV expansion. Thus we start with 
the analysis of the two-point function S(x, x; s) in weakly curved spacetime. 

For this purpose it is convenient to introduce the Riemann normal coordinate system which is a 
coordinate system with the affine connection Y a ^ which vanishes at least locally (i.e., locally inertial 
frame). At the origin xq of the coordinates the metric tensor g M „ and the affine connection T a ^ have 
the property 

rvM = 0, (58) 
where ri^ v is diag(l, — 1, — 1, —1). Near the origin xq the metric tensor is expanded to be p5| 

g„u{y) = r,^ + lRl a ue(y-xor(y-x o ) + O(R.^R 2 ), 

9(y) = -l~\R%{y~x ) a {y-x Y + 0{R.^R 2 ), (59) 

where index for -R° a „^ and RP a p designates the tensor at the origin xq. 

By the use of the Riemann normal coordinate expansion ( |59| ) we calculate the spinor two-point 
function S(x, X] s) in the approximation of slowly varying curvature where we neglect any terms involving 
derivatives of the metric higher than second order. According to the method developed by Parker and 
Tomsp6| the two-point function S(x, x; s) is expanded asymptotically around R° = 0. We introduce the 
bispinor function G(xo,y; s) defined by 

(z 7ai V^ + s)G(x , y; s) = S(x ,y; s) . (60) 

It satisfies the following equation : 

(-V^ - | - s 2 ) G{x , y; s) = j=S D (x , y) . (61) 

7 Since the assumption a c (x) = a (constant) in (S) is not always acceptable in general curved spacetimes, it is necessary 
to evaluate the effective action instead of the effective potential to study the phase structure. Here we restrict ourselves 
to the static homogeneous spacetime where the assumption is acceptable. 
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Using the Riemann normal coordinate expansion ( p9| ) we expand the first term on the left-hand side 
of Eq. @ and find 



rTd^du + \R\p(y - x ) a (y - a*) V^/A 



-|fl OM a(l/-aro) a 0p 



(62) 



1 



G(x ,y;s) , 



where 



(63) 



and Latin indices a and 6 are vierbein indices. Here we keep only terms independent of the curvature 
R and terms linear in R. Inserting the Eq. ( p2| ) into Eq. ( |6l| ) and performing the following Fourier 
transformation, 



G(x ,y;s) = 
we find that Eq. ( |6l| ) reduces to 



(27T) D ' 



-jp(xo-y) 



y; s) , 



UI/ 1 nil 2 ^ n/i 9 1 ntt/3 9 9 

^ - 4* ~ S ~ 3 R ' 6 R 71 P ^W a W, 



1 



i R a 0P » Pv dp a dpp 1 6 



d d 1 a/3 2 d d 
R s 



1 dO^ ^ 



G(p,y;s) = 1 . 



(64) 



(65) 



From Eq. ( pty G(p, y) is found to be 

1 1 BP 



G(p, y; s) = 



. 2 R^ v p,p u 

p 2 - s 2 12 {p 2 - s 2 ) 2 3 (j> 2 - s 2 ) 3 



+ 0(i?° (i? ) 2 ). 



(66) 



Inserting Eqs. ( p4| ) and (|66j ) into Eq. (|60| ) the two-point function S(xo,y;s) in the weak curvature 
expansion is obtained: 0] 



S(ar , y;s) = (^V^ + s) 



d D p 

(27T) D< 



-ip(x -j/) 



G(p,y) 



" ^ ptp(xo-y) 

(2tt) d 



7 Q p a + s J_^ 7 Q p a + s 



p2 _ s 2 



12 (p 2 



(p 2 



T7 17 cda^P 



1 



(p 2 — S 2 ) 2 



+0(i^,(i?T) 



(67) 



It may be easily checked that in the limit R° — > the spinor two-point function in Minkowski space (|30| 
is reproduced. 
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3.2 Effective potential in weakly curved spacetime 

We would like to discuss the curvature induced phase transition in weakly curved spacetime. For this 
purpose we evaluate the effective potential. In the Gross-Neveu type model the effective potential is 
given by Eq. (|25| ) in the leading order of the 1/N expansion. Substituting Eq. (|67j ) to Eq. (|2^) we 
obtain the effective potential in weakly curved spacetime up to terms linear in R 



V(a) = -La 2 



d D p 



ztrl / ds-. 

/o (2tt) d 



_£ Lr 1 

- s 2 12 (p 2 - 



s 2 f 



s 2 f 



(68) 



Here we neglect the index of R° for simplicity. In order to regularize the divergent integral in the 
second term of the right-hand side of Eq. (|6|) we apply the dimensional regularization. Integrating over 
p and s in arbitrary spacetime dimensions we obtain the effective potential 



V(a) = V (a) + V R (a) + 0{R,^R 2 ) . 



(69) 



where Vb(cr) is the effective potential at R 
linear in || 

V R {a) 



given in Eq. 



trl R 



r i 



and V R (a) is the effective potential 



(70) 



(4tt) d / 2 24 

As is mentioned in §2.3 the effective potential Vb(cr) is divergent in two and four dimensions. We 
apply the same renormalization condition as the one in Eq. (|3^) and obtain the renormalized effective 
potential as in the R — case. Replacing the bare coupling constant Ao with the renormalized one A in 
Eq. ( p3| ) we obtain the renormalized effective potential, 



V(a) 



2A 



~ 2 ,,°- 2 



trl 



trl . ( D 



-r 1 1 - — 

{A-k) d ' 2 D \ 2 



trl R 

(4tt) d / 2 24 



J 2 ,, n ~ 2 



D 

1 

2 



(71) 



Before the discussion of the phase structure we consider the two-, three- and four-dimensional limits 
of the effective potential (Fll). Taking the two dimensional limit, D — > 2, we get 



V u=2 {a) 



1 

2A 



trl 
"8V 



-3 + In 



trl R 



D 



7 + In 4-7T 



In I 2 



(72) 



It is different from the expression obtained in Ref. |pl|-F| The effective potential V R (a) (i.e. the third 
term on the right-hand side in the Eq. 

8 Differentiating Eq. ([r2|) we find 



2J)) is also divergent in two dimensions. For D — 2 the 



dV l 



On the other hand Eq. (4.3) in Ref. reads 

dV D = 2 



trl 

47T 



trl 

4-7T 



- In 



In I — 



12a 2 
R 

' 4(T 2 



(73) 



(74) 



in the weak curvature limit where we adopt the notation in the present paper. The common renormalization condition is 
adopted in the above two cases. As is mentioned in §2 the four-fermion theory is renormalizable. The theory is independent 
of the regularization method. Thus the above two results should agree with each other. The solution of the equation for 
the evolution operator in Ref. p]| seems not to be the general solution but a special solution. 
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divergence comes from the lower momentum region of the p-integral in Eq. ( |68| ) . To show it we cut off 
the higher and lower momentum regions of the divergent integral in Eq. (pq) and find 



1 



2A 

trli? 

" 48tt 

-O 



trl 
"8V 



1 



■In 



1 + In 

Atr 
a 



A 



uv 



A 2 



A 2 



A 2 



A 2 



A 2 ' 

ly UV 



(75) 



where Auv is the ultraviolet cut-off scale and Am is the infrared cut-off scale of the divergent integral. 
The ultraviolet divergence is cancelled out after renormalization and the terms independent of R reduce 
to Eq. (|3^). As is seen in Eq. ( f75|) the infrared divergence appears in the terms linear in R for D = 2. 
Because of the infrared divergence the symmetry restoration always occurs for any positive values of the 
spacetime curvature R whereas only a broken phase is realized for negative R. 
Taking the three-dimensional limit D — > 3 we find p3[ 



V u=i (o) 
M 3 
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2A 



trl 

47T 



3 U 



trl R a 

96ir /i 2 /i 



(76) 



If we take the four-dimensional limit D — » 4, the effective potential (|7l| ) reduces to 
^= 4 (q) 

r -i /,\ 2 i /„\ 4S i 

(77) 

V/V J \M/ J 

trl 1 R 

+ 4(47) 2 6^ 2 

where C^ v is defined in Eq. (|38|). 

To see the relationship between C^j v and a cut-off parameter we calculate the effective potential for 
D = 4 by using the cut-off regularization. Performing the integration in Eq. ( p8| ) with the cut-off A in 
the higher momentum region one gets fl22| 




V 



(a) = V»=\a) + 



trl R 

(4tt) 2 6" 



hi 
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- 1 



O 



A 2 



(78) 



where Vb(cr) is given by Eq. (|4C|). Applying the renormalization condition given by Eq. fl32j ) Vo(cr) 
reduces t o (^2] ). After the renormalization we find that there is a correspondence between two results 
(77) and (|78|), if we make a replacement pj 



a 



R 



1 



div 



6 6/i 2 — <t 



i n ; 



A 2 



(79) 



We recognize that the previous correspondence rule ( [43] ) is now modified by the term including curvature 
inEq. ©. 

Now we have succeeded to expand the effective potential asymptotically in terms of the spacetime 
curvature R by using the Riemann normal coordinate. Starting from the effective potential ( |7l| ) we 
investigate the phase structure of the theory in weakly curved spacetime in the following subsection. 
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3.3 Curvature induced phase transition 



The vacuum expectation value is determined by observing the minimum of the effective potential. In 
weakly curved spacetime the effective potential is modified as is mentioned above. Thus the vacuum 
expectation value of a will be changed by the curvature effect. We expect that the phase transition takes 
place by varying the spacetime curvature R. 



3.3.1 Phase structure for A > Acr 

In Minkowski spacetime, as is shown in §2.3, the system is in broken phase for A > Acr- Here in the 
present subsection we fix the coupling constant A larger than the critical value Acr and see whether the 
chiral symmetry is restored by the curvature effect. To study the phase structure in curved spacetime 
we evaluate the effective potential ( |7l| ) numerically with the varying curvature. In Fig. 6 we present the 
typical behavior of the effective potential ( |7l"| ) for several values of the curvature in the case D — 2.5 and 
D = 3.5. Here we adopted the formula trl = 2 D/2 . We find that the chiral symmetry is restored as R is 
increased with A fixed. As can be seen in Fig. 6 the phase transition induced by curvature effects is of 
the first order. On the other hand no phase transition takes place and the vacuum expectation value of 
the auxiliary field simply increases for the negative curvature as its absolute value increases. We observe 
the similar behavior of the effective potential in the spacetime dimensions, 2 < D < 4. 

In the leading order of the 1/N expansion the dynamical fermion mass is equal to the vacuum 
expectation value of the auxiliary field (<j) . The dynamical fermion mass is given through the minimum 
of the effective potential. The necessary condition for the minimum is given by the gap equation : 



dV(a) 



Oct 



= . (80) 



The non-trivial solution of the gap equation corresponds to the dynamical mass of the fermion. Inserting 
the Eqs. ( fn] ) and (^) into Eq. (^) the non-trivial solution of the gap equation is expressed as 

m °->-m°-* + ^(l-£.) m °-* = 0. (81) 



12 V 2, 

In Fig. 7 we plot the non-trivial solution of the gap equation as a function of the spacetime dimension 
D. 

Since the phase transition induced by curvature effects is of the first-order, two independent solutions 
appear for a small positive curvature. The larger solution corresponds to the local minimum and the 
smaller solution represents the first extremum of the effective potential. 

At the critical point the two local minima of the effective potential acquire the same value which is 
zero. Thus the critical point is obtained by the solution of Eq. (|8l]) supplemented by the condition 

V(m) =V(0) =0. (82) 



We solve Eqs. (J8l|) and (p2[) and obtain the critical curvature,) 24 



R cr =Q(D-2)[ ^- D) ) ml, (83) 



and the mass gap at the critical point, |24| 

{D(4-D)\ 1/(d - 2} , . 

TOcr= (^~l ) m °' ( ) 

If the spacetime curvature R is less than the critical value Rcr, the vacuum expectation value of the 
auxiliary field (a) is given by the larger solution of Eq. (|8l|) . In the case that the spacetime curvature 
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Figure 6: Behavior of the effective potential is shown at D = 2.5 and D — 3.5 for fixed A (> A C r) with 
the varying curvature where R C r = 6(D - 2)(£>(4 - D)/4)( 4 " D) / (D - 2 'm^ > 0. 
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Figure 9: Critical curvature as a function of dimension D in weakly curved spacetime. 



R is larger than the critical value i?cr, the chiral symmetry is restored and the vacuum expectation 
value of the auxiliary field (a) disappears. The behavior of the dynamical fermion mass as a function 
of R is given in Fig. 8. We observe that the mass gap appears at R — i? C r and disappears in the 
four-dimensional limit. For some special values of D Eq. (p^) simplifies: 

R cr = ; D = 2 , 

3, (85) 
4. 

In Fig. 9 we show the critical curvature i?cr as a function of the spacetime dimensions D. In two 
dimensions the weak curvature expansion shows that the broken chiral symmetry is restored by the 
effect of an infrared divergence for any positive curvature. The critical curvature i?cr =2 is thus zero. 

3.3.2 Phase structure for A < Acr 

For A < Acr the system is in symmetric phase in Minkowski spacetime. In the present subsection we in- 
vestigate whether the chiral symmetry is broken down by the curvature effect for A < Acr- We introduce, 
for convenience, the scale m' defined in Eq. (|44| ) and calculate the effective potential numerically. In Fig. 
10 the behavior of the effective potential (|7l|) is presented with the varying curvature at D = 2.5 and 
D = 3.5. We observe the second order phase transition as the curvature decreases. The chiral symmetry 
is always broken down for the negative curvature even if the coupling constant A is very small. Thus 
the critical curvature is given by i?cr = in the whole range of the spacetime dimensions D considered 
here, 2 < D < 4. 

In the case with A < Acr Eq. (|8l| ) is modified to be 

- m' °- 2 - m D - 2 + = . (86) 



-Rcr = ! D = 
Rci = l2ml; D = 
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Figure 10: The behavior of the effective potential is shown at D = 2.5 and D = 3.5 for fixed A (< A cr ) 
with the varying curvature where K = 6(D - 2)(D(4 - D) / A)^- D)nD - 2) m' 2 > 0. 
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Figure 11: Solutions of the gap equation for fixed A smaller than Acr &t D = 2.5, 3, 3.5. 



The dynamical mass of the fermion is given by the solution of the Eq. (86). In Fig. 11 we present 
the dynamical mass of the fermion as a function of the spacetime dimensions D. As is expected for the 
second order phase transition, the dynamical mass of the fermion smoothly disappears at R — Rcr = 0. 



3.3.3 Phase structure in four dimensions 

We are interested in applying the results to the critical phenomena at the early stage of the universe. 
Thus we focus on the study of the phase structure in four dimensions here. 

In four dimensions the four-fermion models are not renormalizable. In Minkowski spacetime the 
correspondence ( fl3] ) does not depend on a . Thus both the dimensional and cut-off regularization describe 
the equivalent theory. The differences are removed by the finite renormalization. Therefore it is not 
necessary to investigate the theories defined by the different regularizations in Minkowski spacetime. 

In curved spacetime the correspondence ( f79| ) depends on a and the curvature R. Thus above results 
in weakly curved spacetime may depend on the regularization methods. We regard the theory for 
D = 4 — e with e sufficiently small positive as a regularization of the one in four dimensions. In the 
previous subsections we found that the broken chiral symmetry was restored for a sufficiently large 
positive curvature R > i?cr and the phase transition induced by curvature effects was of the first order 
for a positive finite e. 

To see whether the phase structure depends on the regularization method we define the theory by the 
cut-off regularization. [^|| In this case the effective potential is given by Eqs. ( f40| ) and (j7§) and the gap 
equation defined by Eq. (|o|) reads[] 
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trl 
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(87) 



The dynamical fermion mass is obtained by analyzing the gap equation (|87|). In Fig. 12 the dynamical 
mass of the fermion is plotted as a function of the curvature R for fixed A . The coupling constant 

9 In the present subsection we do not adopt the renormalization procedure. It is not necessary in four dimensions because 
the theory is not renormalizable. 
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Ao is kept in the range Ao > Acr = in 2 /A 2 . The behavior of the dynamical mass as shown in Fig. 
12 is characteristic for the relatively small coupling constant, Acr < Ao < 2Acr- For larger coupling, 
Ao > 2Acr, the behavior near the critical point R = i?cr is quite different from the one in Fig. 12: the 
curve representing the dynamical mass is bent upward near R = i?cr- At any rate there is observed the 
gap in the dynamical mass at the critical point R = Rqt reflecting the nature of the first order phase 
transition. By the direct numerical analysis we find that there is no gap at R = Rqt — if Ao < Acr and 
so the phase transition is of the second order. 

It is possible to obtain the critical value of the curvature R and the coupling constant Ao by observing 
the behavior of the effective potential. The critical values i?cr and A^ r constitute a critical curve in the 
R — Ao plane as shown in Fig. 13. It is found from Fig. 13 that for large positive R the chiral symmetry 
is restored even if A is kept in the region of the broken phase for R = 0. On the other hand the chiral 
symmetry is always broken for negative R irrespective of the value of A. Thus the similar phase structure 
was found in the theory defined by the cut-off regularization. 

To summarize this section we have employed the weak curvature expansion to evaluate the effective 
potential in the leading order of the 1 /N expansion. We assumed that the spacetime is weakly curved 
and have investigated the phase structure. 

Starting from the theory with broken chiral symmetry for vanishing R we calculated the effective 
potential for finite R in arbitrary dimensions 2 < D < 4. We found that the chiral symmetry is restored 
for a large positive curvature.^] The critical curvature dividing the symmetric phase and broken phase is 
obtained analytically. Starting from the theory with chiral symmetry for vanishing R it was found that 
the symmetry is broken down in the spacetime with the negative curvature. Therefore only a broken 
phase occurs for the negative i?[^2|, |24| 



4 Exact solutions in special spacetimes 

In the previous section the curvature induced phase transition was discussed by using the weak curvature 
approximation. It is found that the broken chiral symmetry is restored for a large positive curvature. 
But the large positive curvature may spoil the validity of the weak curvature approximation. To study 
the phase transition in the large positive curvature in more confident manner we need to avoid any 
approximation in dealing with the spacetime curvature. In some specific spacetimes the effective po- 
tential is calculable rigorously in the leading order of the 1/N expansion without using weak curvature 
approximation because of the symmetry of the spacetime. 

In the present section we try to find the effective potential of the four-fermion model without making 
any approximation in the spacetime curvature. For this purpose we restrict ourselves to specific space- 
times, i.e., maximally symmetric spacetime (S D , H D ) |27| [2^, ||^, |3(| and Einstein universe (R® S ^ 1 ) 
pi] , p2] , p3| , p^l and calculate the effective potential in the leading order of the 1/N expansion. The 
phase structure of the four-fermion model is obtained by using the same method developed in §2 and 
§3. To justify the weak curvature approximation we compare the result with the one obtained in §3. It 
is shown that the weak curvature approximation gives the exact result [p4] |28| , [m] in the limit D — > 4 
or A -> oo. We follow here Refs. @ |§ and @. 

4.1 Spinor two-point function in S D and H D 

As we have seen in §2 the effective potential for the composite field in the leading order of the 1/N 
expansion for the Gross-Neveu type model is described by the two-point function S{x 1 x; s) of the massive 
free fermion in curved spacetime. Thus we start with the analysis of the two-point function S(x, y; s) in 
the maximally symmetric spacetime. 

In the maximally symmetric spacetime S D and H D (2 < D < 4) the exact expression of the two-point 
function is known without making any approximation with respect to the spacetime curvature. [[36] |37| , 

10 Large positive curvature may spoil the validity of the weak curvature expansion. This point will be discussed in §4. 
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p8| . ^9) We then restrict ourselves to the maximally symmetric spacetime and closely follow the method 
developed by Camporesi.]39| Here we consider the manifolds S D and H D . The manifold S D is defined 
by the metric 

ds 2 = a 2 {d6 2 + sin 2 6dn o _ 1 ) , (88) 
where dVl D _ l is the metric on a unit sphere S ^ 1 while the manifold H D is defined by 

ds 2 = a 2 (d8 2 + sinh 2 edVt D ^) . (89) 

The manifold S D and H D are constant curvature spacetimes with positive and negative curvature 

R = ±D(D - l)a~ 2 , (90) 

respectively (2 < D < 4). 

The spinor two-point function S(x,y;s) is defined by the Dirac equation 

(f + s)S(x,y;s) = -^-S D (x,y), (91) 
V9 

where S D (x,y) is the Dirac delta function in the maximally symmetric spacetime. We introduce the 
bispinor function G defined by 

(f- S )G(x,y;s)=S(x,y;s). (92) 
According to Eq. ( |9l| ) G(x, y; s) satisfies the following equation, 

(ff- s 2 )G(x,y;s) = -^-S D (x,y). (93) 
On S D and H D we rewrite Eq. ( p3| ) in the following form 

U D - I - ,s 2 j G(x, y- s) = -^=* D & V) - (94) 

where O d is the Laplacian on the maximally symmetric spacetime. 
The general form of the Green function G(x, y; s) is written asfjof 

G(x,y;s)=U(x,y)g D (l), (95) 

where U is a matrix in the spinor indices, g D is a scalar function depending only on I, I = a6 which is the 
geodesic distance between x and y on the maximally symmetric spacetime, rii is a unit vector tangent 
to the geodesic rii = VJ. Inserting Eq. (|95|) into Eq. (|94|) we get 



0, (96) 

where we restrict ourselves to the region I ^ 0. To evaluate Eq. (|96|) we have to calculate the covariant 
derivative of U and n^. U is the operator which makes parallel transport of the spinor at point x along 
the geodesic to point y. Thus the operator U must satisfy the following parallel transport equations: ||37| 

n l ViU = 0, 

(97) 

U(x,x) = 1. 
To evaluate the second derivative of U we set 

ViU = ViU. (98) 



29 



From the integrability condition |3(| on Vi, 

ViVj-VjVi-[Vi,Vi] 



R 



D(D-l) 13 ' 

and the parallel transport equation (f)7\) we easily find that 



1 



Vi = tan — (Jijn 3 ; on S , 



2a 



Vi = — tanh 
a 



<7ijn 3 : on H 



D 



(99) 



(100) 



(101) 



where ov,- are the antisymmetric tensors constructed by the Dirac gamma matrices, = i[7i,7j]- 
To find Vi we have used the fact that the maximally symmetric bitensors are represented as a sum of 
products of rii and g^ with coefficients which are functions only of L ps} After some calculations we get 
the Laplacian acting on U 



u D U = - 



D- 1 

D - I 

4a 2 



tan 



■ tanh" 



I 

2a 

I 

2a 



U, onS D , 
U, onH D . 



(102) 



The derivative of rii is also the maximally symmetric bitcnsor |pq| and found to be 



ViUj = - cot ( - ) (//,, - //;/),). 



S D , 



1 



(103) 



ViUj = — coth — I (gij — riiUj) , on H 



Therefore Eq. ( Pq) reads 



D - 1 



cot - } di - 



D- 1 

4a 2 



2(i 



R 



tan^ — - — - s z ) g D = , on S 



dr 



D- 1 



■ coth 



di- 



D - 1 



4a 2 



■ tanh 



2 a 



R 

- ~i - s z ) 9n = , on H 1 



(104) 



(105) 



We define the functions Hsd{1) and Hhd{1) by g D (l) — cos(l/2a) Hsd{1) an d g D (l) = cosh (I /2a) ha d (I) 
respectively and make a change of variable by z = cos 2 {I /2a) in Eq. (104) and z' — cosh 2 (I /2a) in Eq. 
(105). We then find that Eqs. (104) and ( |105| ) are rewritten in the forms of hypergeometric differential 
equations: 



z{\ - z)d z 2 + 



z\\-z')d z> 2 + 



D + 2 



D + 2 



(D + l)z d. 



2 2 

s a 



Iisd(z) = 0. 



- (D + l)z' d z . - 



D 2 



2 2 

+ s a 



h H n(z') = 0. 



(106) 



(107) 



Noting that the Green functions are regular at the point / = an and fall off for I — > oo we write the 

(108) 



solutions of Eqs. (106) and (107) by the hypergeometric function 

D 



hsD(z) = c S dF 



D D + 2 

— h tsa, isa, ; z 

2 2 2 



h HD (z r ) = c HD {-z')- D,2 - sa F + sa, sa, 2sa + 1; - f 



(109) 
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As we remained in the region where I ^ the normalization constants csd and chd are yet undetermined. 
To obtain csd and crd we consider the singularity of G(x, y; s) in the limit I — > 0, 



r 



G 



G 



csd ~ 



D + 2 



D-2 



D 



r h isa r isa 



D 



CffD(-l) 



_ U ~D/2- 



T(2sa+ 1)T 



I 

2a 
D-2 



D 



+ sa ) r (sa) 



2-D 



(110) 



and compare them with the singularity of the Green function in the flat spacetime. This procedure is 
justified because the singularity on a curved spacetime background has the same structure as that in 
the flat spacetime. For I ~ the Green function in the flat spacetime behaves as 



G flat (0 



1 



4tt £ >/ 2 



D-2 



I 2 



(111) 



Comparing Eq. flllO]) with Eq. (Ill) the over-all factors csd and chd are obtained: 



,2-n 



CSD 



CHD 



D 



isa r 



D 



{An) 1 



(_-i\D/2+sa_^ 

' ' (4tt) d 



D + 2 



j ; , r ( — +sa) Visa) 



T(2sa+l) 

Inserting the Eqs. (|5g), (pM and (pT2|) into Eq. @ we find on S D 



(112) 



G(x,y;s) = U{x,y) 



D 



F h isa r isa 



D 



(4tt) d / 



X cos 



F 



D 



D + 2 



D 



(113) 



+ isa. isa, ■ 

2 '2 



D + 2 



: cos 



and on H 



G(x,y;s) 



U(x,y) 



a*-° F (f + * a ) r M 
( 47r p/2 T(2sa+1) 

n l-D-2sa 



(114) 



cosh 



F 



D 



sa. sa, 2sa + 1; cosh 



2rt 



Thus the Green functions G(x, y; s) on the maximally symmetric spacetime are obtained. 

The spinor two-point function S(x,y;s) is derived from the Green function G{x,y;s). From the Eq. 
© we get 



5 = (7 l V. ( - s)Ug D 

- , ii 7 j d t - D 1 tan ( ^- ) ) g D - sll g D , on S D , 



jirfU di + 



D — 1 
2a 



tanh ( — 
2a 



— sUg D , on H 



D 



(115) 
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Substituting Eqs. (113) and (114) in Eq. (|115|) the spinor two-point function S(x,y;s) is obtained 



S(x,y;s) 



on S D and 



S(x,y;s) 



D 



r — + isa r — — isa 



D 



(4tt) d / 



sU (x, y) cos 



D 



D + 2 



2a 



+j i n'U(x,y) — sin 
la 



F 



2a 



D 



D D + 2 2 
isa, isa, ; cos 



D 



D 



D 



isa, isa, — ; cos 

2 2 



!_ D r(^ + sa )r( sa + i) 



(471-)" 



/2 



T(2sa+ 1) 



cosh 



l-D-2sa 



2a 



D 



[/(a;, y) cosh I — ] F I — + sa, sa, 2sa + 1; cosh I — 



2a 



+-f i n 1 U(x,y)smh ( — 



D 



• ./• j ™ + >•(/. «» ■+ .1 . 2s/i ■+ .1 : <•(»!) - ( J- 



(116) 



(117) 



on H D . According to the anticommutation relation of spinor helds the two-point function (116) satisfies 
the antiperiodic boundary condition S(l) = —S(l + 2i:na) where n is an arbitrary integer. 

We succeeded to calculate the spinor two-point functions without making any approximation in the 
spacetime curvature. Using these functions we evaluate the exact expression of the effective potential in 
curved spacetime in the following subsections. 



4.2 de Sitter background (S D ) 

First we consider the four-fermion model in de Sitter space [^7[ |28|, ^9). The D-dimensional de Sitter 
space is represented as a hyperboloid, 



t2 t 2 
?0 - SI 



f 2 



(118) 



embedded in the (D + l)-dimensional Minkowski space. It is one of the maximally symmetric spacetime. 
The de Sitter space is a constant curvature spacetime with curvature 



R = D(D - l)a~ 



(119) 



We consider the manifold S D as a Euclidean analog of the Z?-dimensional de Sitter space. As is 
mentioned in the previous subsection the spinor two-point function S(x,y; s) is known on the manifold 
S D . trS(x, x; s) is required in evaluating the effective potential: 



trS(x,x;s) = - (4 ^ o/2 T| I -- 



D 



r — + isa r — — isa 



D 



r (1 + isa) r (1 — isa) 



(120) 
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Performing the Wick rotation and inserting Eq. (120) into Eq. 
potential is obtained. p7J 



J) our final expression of the effective 



V(a) 



trla 2 



2An 



(47T) 



D/2 



-r i 



r 



sds- 



D 
~2 



D 



T{\ + isa)T{\ - isa) 



(121) 



Equation (121) is an exact expression of the effective potential for the model of four-fermion interactions 
in de Sitter space in the leading order of the 1/N expansion. Substituting the renormalized coupling 
constant A defined by Eq. (|3^) in the Eq. (121) we find the renormalized expression of the effective 
potential V(a) in de Sitter space, 



V(a) 



1 

2A 



a fj, 



trl 



2(4tt) d 



/2 



(D-i)r l- 



D 



trio 2 



(47T)D/a 

Note that Eq. jl2^ ) reduces to 

y D = 2 (q) 

M 2 



r l- 



J 1 ,, n ~ 2 



1 isa 

2 



T (1 + isa) r (1 — isa) 



(122) 



1 

2A ~ 
trl 

47r,u 2 



S ( ln(Q/i)2 



ciss ["0(1 



2) , 
J VM, 

isa) + "0(1 



(123) 



in two dimensions. Expanding Eq. (122) asymptotically around 1/a = (weak curvature expansion) 
the effective potential (|34|) is reproduced. In Fig. 14 the behavior of the effective potential given by 
Eq. (122) is illustrated in the case D = 2.5 for several typical values of the curvature. It is observed 
in Fig. 14 that, if A < Acr, the theory is always in the symmetric phase as the curvature changes 
while, if A > Acr, the symmetry restoration takes place as the curvature exceeds its critical value. This 
observation remains true if the spacetime dimension is arbitrarily changed. 

To discuss the dynamical mass of the fermion we study the minimum of the effective potential more 



precisely. A necessary condition for the minimum is given by the gap equation fl80|). Inserting Eq. (122) 
into Eq. (§3) we find Q 



trl- 



D - 1 



(47r)D/2 

trl(a/i) 2 - D 
(4tt) d 



-r i 



D 
~2 



r i 



D 
~2~ 



D 



+ ima r 



D 



r (1 + ima) r (1 — ima) 



(124) 



If the coupling constant A is no less than a critical value A cr , the gap equation allows a nontrivial solution. 



In Fig. 15 we present the solution of the gap equation (124). The dynamical fermion mass smoothly 
disappears as the curvature increases. 



Taking the two-dimensional limit Eq. (124) reduces to 

1 trl 

A ~ 4vr 

For three dimensions it reads 



[ln(a 2 ^( 2 ) — ^(1 + ima) — — ima)\ = . 



(125) 



1 

A 



2trl 

47T 



trl 



Aira/i r (1 + ima) T (1 — ima) 



0. 



(126) 



33 





Figure 14: Behavior of the effective potential is shown at D — 2.5 for fixed A with the varying curvature 
where K = 8.23m'o and R CT = 8.23ml 
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Figure 15: Behavior of the dynamical fermion mass in de Sitter space as a function of the curvature R 
at D — 2.0, 2.5, 3.0, 3.5, 4 where too is the dynamical fermion mass in flat spacetime. 



while it reduces to 



1 



3trl 



2\ 2trla~ 2 
A (4tt) 2 V div ~ 3 J ~ (450V 



(4tt)V 



[^div + ln(a 2 /i 2 ) — ip(l + ima) — -0(1 ™ ima)\ = , 



(127) 



in the limit of D — ► 4, where the divergent part is expressed by C f\\ v in Eq. (38). Taking the four- 
dimcnsional limit, D — > 4, the divergent parts of the gap equation (127) is exactly equal to that given 
by Eq. (|8l|). Hence the same critical behavior is observed for the dynamically generated fermion mass 
in Figs. 8 and 13 at D = 4. 

As is seen in Figs. 12 and 13 the phase transition is of the second order. Accordingly by setting 
to = in the gap equation (124) we may derive the equation which determines the critical radius acr, 



1 , D-l 
- + trl- r 

A (4tt) d / 2 



D 



D 



D 



r z — r i - — =0 



Substituting Eqs. 



and 



r 1 trl , . 

V 2 J (^tt) d ' 2 \ 2 J \ 2 

in Eq. (|128|) the critical radius acr is found to be (27 



acr 



1 

to 



2/(D-2) 



(128) 



(129) 



In Fig. 18(a) t he cr itical radius acr is plotted as a function of the spacetime dimension D. By using the 
critical radius ( 129 ) we find the critical curvature: 



i?cr = D(D - l)a cr " 



(130) 
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For some special values of D Eq. ( |130| ) simplifies to: 

R CT = 2e 2 -<m 2 ;D = 2. 



Rcr = ^ml ;D = 3, (131) 



Rc T = 12ml ;D = A, 

where 7 is the Euler constant. In Fig. 18(b) the critical curvature i?cr is presented as a function of the 
spacetime dimension D. 

4.3 Anti-de Sitter background (H D ) 

Next we consider the anti-de Sitter space. It is also the maximally symmetric spacetime with negative 
curvature 

R =-D(D- l)a~ 2 . (132) 

We consider the manifold H D as a Euclidean analog of the D-dimensional anti-de Sitter space. 

As is discussed in §4.1 the spinor two-point function S(x, x;s) is known in H D and trS(x,x;s) is 
given by 

'd \ „ / zr 



trla^U +H r ( 1 - 2 

: 1 - s ' ( ' • ; ' : * s = - i&m f d ^ • 

r [ 1 - — + sa 



Performing the Wick rotation and inserting Eq. (132) into Eq. ( J25| ) l ne effective potential of the 
four-fermion model reads 

Thus we obtain the exact effective potential ( |l34| ) in anti-de Sitter space in the leading order of the 1/N 
expansion. Applying the same renormalization condition (32) in Eq. (134) we find the renormalized 
expression of the effective potential V(a) in anti-de Sitter space, 

Via) = — f rV' 2 + (D - l)r ( 1 - - VV~ 3 

V ; 2A P 2(4tt) d / 2V ' \ 2 J P 

Expanding Eq. (|135|) asymptotically about 1/a = (weak curvature expansion) the effective potential 
(34h is reproduced. To study the phase structure in anti-de Sitter space we evaluate the effective potential 

In order to illustrate phase structure of four-fermion model in anti-de Sitter space we will consider 
two-dimensional case only below. That will be enough for our purposes because as we will see negative 
curvature always supports chiral symmetry breaking. Hence, there is no need to make the exposition of 
cases for all D as it was above. Note that the fact that negative curvature supports chiral symmetry 
breaking and only broken phase survives has been already mentioned few times when we have discussed 
weak curvature expansion. 
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Figure 16: Behavior of the effective potential in H 2 



In two dimensions the effective potential (135) reduces to ]2 
V D=2 (a) 



2A 47r \ a/i 
trl 



eLss [tf>(l + sa) + V'('Sa)] 



(136) 



Evaluating the effective potential (136) numerically we easily see that only the broken phase is realized 



in two dimensions. In Fig. 16 the behavior of the effective potential given by Eq. ( Jl 36| ) is illustrated. 
Differentiating (|136[) with respect to a we get 



dvD=2 f 1 trl u . 2 2 , , . . . Nl 

d(j = o | - - — [hi(aV) - ^(1 + ^a) - i>{aa)] 



(137) 



A careful study of dV/da shows that due to the presence of the last term in (|137j), a = is never 
stationary for any value of A and finite a. Owing to the fact that 



dV 



D=2 



Oct 



trl 

Ana 



(138) 



the chiral symmetry is always broken in two-dimensional anti-de Sitter space. ||29| 
In arbitrary dimensions (2 < D < 4) we easily find 



dV 



a=0 



trl , n f D\ 

- — : a}- D Y — < 0. 

4tt) d / 2 V 2 



(139) 



Thus the chiral symmetry is always broken down in anti-de Sitter space (i.e., a negative curvature 
spacetime) in arbitrary dimensions 2 < D < 4 as is mentioned in §3. 
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4.4 Einstein universe background (R®S D x ) 

Another universe we can calculate the effective potential without making any approximation in the space- 
time curvature is the static Einstein universe. |34| The D-dimensional Einstein universe is represented 
by the metric 

ds 2 = dr 2 - a 2 (d9 2 + sin 2 6dn D _ 2 ) . (140) 
It is a constant curvature spacetime with curvature 

R=(D-l)(D-2)±. (141) 

Here we consider the manifold R(g> S D ~ 1 as a Euclidean analog of the D-dimensional Einstein universe. 
In two spacetime dimensions a Euclidean analog of the Einstein universe is represented by a cylinder 
with radius a. 

According to a similar analysis developed in §4.1 we are able to solve the Dirac equation in Einstein 
universe. |^5) On S 0-1 we rewrite Eq. (|93| ) in the following form 

(do) 2 + D D _i - - - /) G(x, y; s) = ~^ D {^ v) > (142) 

where Dd-i is the Laplacian on S 13-1 . Performing the Fourier transformation 

G(x, y;s) = J ^e^-^G^, y - x) , (143) 



we rewrite Eq. (142) in the form 



-'/>-; ---(s*+^))G(u,y-x) = -—5 u -"{y). (144) 



Equation (144) is of the same form as the one for the spinor Green function with mass \/ s 2 + ui 2 on 
$d-i^ ^hus explicit expression for the two-point function in i?£g> S 0-1 is given by 

, q/ s trW-° /3-D 

txS(x,x;s) = -» (47r)(o _ 1)/2 T 



djj 



2?r T(l + ia)T(l-ia) 



(145) 



where the parameter a is defined by 

a = a^s 2 +uj 2 . (146) 



Inserting the two-point function (145) into Eq. (|25| ) we obtain the effective potential for the four- 
fermion model in Einstein universe in the leading order of the 1/N expansion. |]34[ 

1 2 trlo3-° /3-D 
V{<7) = 2A^ (^-^ T [— 

1 — h 1a 1 — ia 



srfs / — v " ' V " ^. (147) 

' 2tt r(l+to)r(l-ia) 



38 



We apply the renormalization condition ( |32| ) and obtain the renormalized effective potential by replacing 
the coupling constant Aq with the renormalized one A defined by Eq. 



V{a) = 



~2,,D-2 

2A* " 



trla J 
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2(4tt) d / 
"i-D 
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(148) 



Expanding Eq. (|148|) asymptotically about 1/a = (weak curvature expansion) the effective potential 
(34) is reproduced again. For convenience in a numerical calculation we rewrite Eq. (148) as 
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(4tt)( d 
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2vr 



r(l+fa)r(l - ia) 



By the use of the expression of the effective potential (149) we may study the behavior of the effective 
potential as a function of a through the numerical integrations. In the numerical integration in to we 
introduced a suitable upper and lower bound and performed the numerical integration between these 
bounds. By assuming the sufficiently large absolute value of these bounds and checking the stability of 
the integral under the change of the bounds we obtained the numerical value for V(a) for each value of 
a with A and a kept fixed. 

A behavior of the effective potential is similar to that in de Sitter space illustrated in Fig. 14. If 
A < A cr , the theory is always in the symmetric phase as the curvature changes while, if A > Acr ; the 
symmetry restoration takes place as the curvature exceeds its critical value. Only the second order phase 
transition occurs as the curvature exceeds its critical value. |f34 

The dynamical fermion mass is given by the gap equation (8(j): 



V- + — r ( (^L\[^ + (d- 2) M 2 ](^ 2 + M 2 )(— ^ 

y (4tt)( d - 1 )/ 2 V 2 J J 27rV V ^ U ^' 



rl D-l \ r D-l . ^ 
1 h ia 1 ia 



3-D 



r(l + ia)T(l - ia) 



= 0. 



(150) 



Figure 17 represents the behavior of the dynamical fermion mass which is obtained by solving Eq. 
(150) numerically for D = 2.0,2,5,3.0 and 3.5. In two dimensions Einstein universe becomes a flat 



spacetimc, R = 0. To see the mass dependence of the radius a for D — 2 we draw the dynamical fermion 
mass as a function of (am )~ 2 instead of the curvature R. 

T he cr itical radius is given by the massless limit of the gap equation. Taking the massless limit in 
Eq. ( |L50[) , we can perform the integration over u> analytically and find 
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(151) 
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Figure 17: Behavior of the dynamical fermion mass in Einstein universe as a function of (amo) 2 at 
D = 2.0, 2.5, 3.0, 3.5 where mo is the dynamical fermion mass in flat spacetime. 



Taking into account Eq. 



i]) with Eq. (J49[) we rewrite Eq. (151) in the following form 

l/CD-2) 



acr 
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m 



V 2 



(152) 



Solving Eq. (152) the critical radius is obtained in an arbitrary dimension. In Fig. 18(a) we show the 
dependence of the critical radius to dimension D. It should be noted that Eq. ( |l52| ) reduces to 



acr 



2toq ' 



(153) 



i n tw o dimensions which reproduces the result obtained in Refs. [J40| and Using the relationship 
(141), the critical curvature is obtained from the critical radius. For some special values of D the critical 



curvature i?cr simplifies to 



^cr 
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;D 


-2, 
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= 8ml 
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-3, 
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= 12m 2 , 


;D 


= 4. 



(154) 



The critical curvature i?cr is plotted in Fig. 18(b) as a function of the spacetime dimension D. It is 
clearly seen in Fig. 18(b) that three lines of the critical curvature reach the same value 12m 2 , at D — > 4. 
For D = 2 the critical curvature obtained by weakly curvature approximation is exactly equal to that in 
Einstein universe. As is explained in the last section the chiral symmetry is restored by the effect of an 
infra-red divergence for any finite curvature in two-dimensional weakly curved spacetime. The critical 
curvature is thus zero. In two-dimensional Einstein universe it is obtained that the critical curvature 
^?cr = 0. The situation is, however, different from that in weakly curved spacetime. By definition two- 
dimensional Einstein universe is a flat spacetime, R = 0. The symmetry restoration is induced by finite 
size effects of the compact and closed space. The spacetime curvature R is not suitable to represent the 
phase structure in two-dimensional Einstein universe. 
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(a) Critical radius a C r as a function of dimension D. 
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(b) Critical curvature i?cr as a function of dimension D. 

Figure 18: Phase diagrams in de Sitter space and Einstein universe. The full and the dashed lines 
represent the exact solution in de Sitter space and Einstein universe respectively. The dotted line 
represents the critical curvature obtained by the weak curvature approximation. 
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Thus we have considered the Gross- Neveu type model as one of the prototype models of the dynamical 
symmetry breaking and investigated the phase transition induced by the curvature effect without making 
any approximation in the spacetime curvature. 

In de Sitter, anti-de Sitter space and Einstein universe the two-point functions have been solved 
exactly. Then the exact expression of the effective potential is obtained in such spacetimes. We calculated 
the renormalized effective potential for finite R in the leading order of the 1 /N expansion in arbitrary 
dimensions 2 < D < 4 . It is found that the broken chiral symmetry is restored at a certain critical 
curvature in de Sitter space and Einstein universe (i.e., positive curvature spacetimes). The phase 
transition from the broken phase to the symmetric phase is of second order. In anti-de Sitter space (i.e., 
negative curvature spacetime) only the broken phase is realized irrespective of the value of A. 

For a negative curvature spacetime the weak curvature approximation gives the exact phase structure, 
Rev — 0, in the whole range of D considered here: 2 < D < 4. The critical curvature i?cr obtained 
by weak curvature approximation is exactly equal to the one obtained in de Sitter space and Einstein 
universe at D = 4 shown in Fig. 18(b). In four dimensions ultraviolet divergences appear in terms 
independent of the curvature R and terms linear in R only. The higher order terms in R of the effective 
potential are ultraviolet finite. Expanding the exact results in de Sitter space and Einstein universe 
asymptotically about R = 0, we obtain the R 2 term of the effective potential in a compact spacetime 
with a constant curvature. 

V(a) = Vo(ct) + V R (a) + V R2 (a) + 0(i? 3 ) , (155) 
where Vo(ct) and V R (a) are given by Eqs. (J3l|) and (|70|) respectively, the R 2 term V R2 (a) reads pj 



V M = -**T (l - °) ^P-2)P-3)(2 + 5J ) 

V ' (4tt) d / 2 V 2 / 5760 D(D-l) V ' 



At the four-dimensional limit the R 2 term (156) reduces to 



V£T\a) _ trl UR 2 / 173 , J 



/i D (4tt) 2 17280 I aiv 66 



^div-^T-ln - • (157) 



The divergent parts in Eq. ( 157 ) appear from the mass singularity at a — > and the normalization 
condition V(0) — 0. Only an infrared divergence appears in the R 2 term ( 157 ). The infrared divergence 



does not appear in the de Sitter space and Einstein universe since the space components are compact 
and closed (i.e. there is an infrared cut-off Am = l/(2a)). The R 2 term has no contribution to the 
symmetry restoration in a compact space (de Sitter space and Einstein universe). Thus the critical 
curvature is determined by the terms involving ultraviolet divergences in four dimensions and the weak 
curvature approximation gives the exact result. Therefore the weak curvature approximation seems to 
be useful near four dimensions. [^| If we use the cut-off regularization in four dimensions we also see the 
similar situation. In the case the weak curvature approximation also gives the exact result at the limit 
A ->■ oo.[|| 

In the other dimensions the critical curvature ( |85| ) is unable to compare with the results obtained in 
de Sitter space and Einstein universe directly, because the global topology of the spacetime may play a 
crucial role for the symmetry restoration. In two-dimensional Einstein universe R<g> S not the curvature 
but the size of the space coordinate a determines the minimum of the effective potential. (See §6.) The 
critical radius acr in Einstein universe is half as many as that in de Sitter space. It seems that the 
number of the compactified directions determine the critical radius. Thus the finite size effect may play 
an essential role in two dimensions P| 

For D = 2 the chiral symmetry may be restored at any value of a in the case of finite JV through the creation of an 
instanton-antiinstanton or kink-antikink condensation as is discussed at finite temperature. [4l| In our method we cannot 
deal with the influence of time or space dependent configurations. Instantons and kinks are, however, suppressed at the 
large N limit and the phase transition can take place in two dimensions. 
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5 Extensions of four-fermion model 



5.1 Gauged NJL model in curved spacetime 

Gauged NJL model maybe considered as one of the most realistic four-fermion models. [Q From one 
side, that is one of the best effective theories to describe QCD. From another side, gauged NJL model 
may well describe SM (standard model) where composite bound states play the role of an elementary 
Higgs fields in the process of dynamical symmetry breaking. 

For a long time the only way to study gauged NJL model has been related with Schwinger-Dyson 
equation. p3|'p^| However in Ref. Q| (following the renormalization group (RG) approach of Ref. ^5|) 
it has been shown that some class of gauge Higgs- Yukawa models in leading order of a modified l/N c 
expansion leads to a well-defined, non-trivial theory. This theory is equivalent to the gauged NJL model 
when the ultraviolet cut-off goes to infinity (using corresponding compositeness conditions). Then, it 
appears the possibility to find the effective potential and to discuss the phase structure of gauged NJL 
model without use of Schwinger-Dyson equations, ft has been done in flat space in Ref. |44j and in 



curved space in Ref. 46 . Below, we will discuss the phase structure of gauged NJL model following 
closely Ref. @. 

5.1.1 Renormalization group equations for the class of gauge Higgs- Yukawa models 

In this subsection we review the class of gauge Higgs- Yukawa models in flat spacetime we are going to 
extend and the approximations we are using in such an analysis. Thereby we closely follow Ref. p4[ 
where, using the RG approach, this model has been considered in detail. 

The Lagrangian of the model in flat spacetime is given by (we use the notations of Ref. Q] ) 

r 1 rta ria^iv , 1 fa _\2 1 ™2 2 ^ 4 

L ra = - j £>,G +2 v m ) - 2 ~ 4 

Nf n; 

+ J2 & i ^ D ^ ^ - y^^i ■ ( 158 ) 

i=l i=l 

Here the gauge group SU(N C ) is chosen, with Nf fermions ipi(i = 1,2, ...,Nf) belonging to the repre- 
sentation R of SU(N C ), a is the scalar field. 

Let us now describe the I /N approximation for the perturbative study of this theory at high energies 
through the RG equations: 

a) The gauge coupling constant is assumed to be small: 

q± « i , (159) 

and the RG equations are considered only in the first non-trivial order on g 2 . 

b) The number of fermions should be large enough: 

N f ~ N c . (160) 

On the same time it is supposed that only nf fermions (jif <C Nf) have large Yukawa couplings, whereas 
the remaining fermions have vanishing Yukawa couplings. 

c) The approach is assumed to be perturbative also in 1/N C , so only the leading order of the 1/N C 
expansion survives; this means that scalar loop contributions should be negligible 



< N c . (161) 



12 It is not known how to formulate consistent Schwinger Dyson equation even in ladder approximation in general curved 
spacetime. 
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Note that above approximation maybe reasonable for the minimal SM. 

Within the above approximation the RG equations for the coupling constants in (5-1) are^| 



dg(t) 



9 3 (t) , 

■ii)- 

dy(t) y(t) . 2 



dt 



(4tt) 2 



dt 



(47T) 



dX(t) _ u y 2 (t) 



dt 



(4tt) 2 



y\t)-cg*{t)] 



m-y 2 (t)) 



(162) 



where b = (11N C - 4T(R)Nf)/3, c = 6 C 2 (R), a = u/4 = 2 njN c . For the fundamental representation 
we have T(R) — 1/2, C2(R) — (N 2 — X)/(2N C ). Here t = ln(/x//xo) and /j,q is the reference scale to discuss 
low energy physics. 

The RG equation for g 2 (respectively a = g 2 /4ir) is solved by 



T)(t) 



g 2 (t) _ a(t) 



2vr 



(163) 



To solve the RG equations for the Yukawa and the scalar couplings it is convenient to introduce the 
following RG invariants: 



With their help we obtain 



h(t) = -rf 
k{t) = -rf 

lf(t) = 
A(t) = 



-l+c/6( t ) 
-l+2c/6 (t) 



1 - 



1 - 



c-b g 2 (t) 
a y 2 (t) 

2c -b X(t) g 2 (t) 



2a y\t) y 2 {t) 



(164) 



g 2 {t) l + M Wb W 

a L 

2a y\t) 



2c -b g 2 (t) 



i + W- 2c ' b {t) 



(165) 



Of course, these solutions are explicitly known only when the values of the RG invariants h and k are 
given. Note that in (5-8) and below we use notations ho and ko for h and k in order to show that they 
are constants (RG invariants). 

As it has been shown in detail in Ref. |Q for the solutions (5-6) ~ (5-8) to be consistent with the 
above approximation (5-2) ~ (5-4), and in order for the theory to be non-trivial one, we should have 



c > b and < ho < oo 



(166) 



That makes the Yukawa coupling to be asymptotically free and the theory to be non-trivial. Further 
analysis shows that the scalar coupling constant is nontrivial in the above approximation only if 



fcn = 



(167) 



Below we will discuss the solutions (5-8) with conditions (5-9) and (5T0) only. 

Let us extend now this model to curved spacetime. It is by now well-known (see Ref. for an 
introduction) that in order to be multiplicatively renormalizable in curved spacetime, the Lagrangian of 
the theory should be: 

L = L m + L ext - hRcr 2 , (168) 



For general discussion of one-loop RG equations in gauge Higgs- Yukawa models, see Refs. 



and 
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where L m is given by (5-1) with a change of flat (partial) derivatives to the corresponding covariant 
derivatives: <9 M — ► Vu, £ is the non-minimal scalar-gravitational coupling constant, and 



-^ext = a-iR 2 + a 2 C 



fiva/3 



a 3 G + a^UR + A - -R . 



(169) 



The Lagrangian £ ex t of the external gravitational field is necessarily introduced in order to have the 
theory to be multiplicatively renormalizable one; a±, ...04, k are gravitational coupling constants. 

Since the RG equations for those coupling constants which are present in flat space do not change in 
curved spacetime,(7| all the discussions of RG equations for g,y and A (as well as the approximations 
introduced above) are also valid here. In addition, the effective coupling constants corresponding to 
£, 01, 04, G appear. However, for our purposes, we only need the RG equation for the coupling 
constant £ which maybe written as following (within the above described approximation), 



dm 

dt 



(4tt) : 



2ay 2 (t) 



Taking into account that 



d „, s d 



, /b{t) m - 1) 



0. 



we may use this RG invariant f(t) to find the solution of (5.13): 

m = l + y 2 (t)v- c/b (t)f , 

where /o ^ is the value of the RG invariant. Note that in the UV-limit (t — > 00) we have: 

,2 



(170) 



(171) 



(172) 



y 2 (t) 



— 9 f v c/b (t) 

a riQ 

fo 



-0, 



(h > 0) 



9 u 

a Iiq 



(173) 



For fo « the scalar-gravitational coupling constant £(£) — > 1/6, i.e., asymptotic conformal invariance 
is realized. f50[ Hence, the RG invariant fo characterizes the deviation from conformal invariance (£ = 
1/6). Different types of UV-behavior of for different GUT's have been listed in Ref. 0; the most 
typical ones are: 



a) £(*)-g , 

b) |f(t)|-^oo 

c) £(i)^£ . 



(174) 



As we can see, unlike the analysis in the case of flat spacetime, there will not appear any restrictions to 
the sign and the value of fo from the study of the UV-asymptotics of £(t). 

Let us discuss now the one-loop effective potential for the theory (5.1) in the approach where we keep 
only terms with the accuracy up to linear curvature, [[i"9[ [5l| i.e., a 2 » \R\. This approach is actually 
equivalent to the one described in §3. Within above approximation we get 



V 



1 
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-fRa 2 
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/4 



- 1 



(175) 



where fip = ya. 

Using such a form for the potential V one may investigate the curvature-induced phase transitions 
which have been discussed in detail in Refs. [[l9[ f?J and |^1[ for different GUT's. 
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5.1.2 Gauged NJL model in curved spacetime 

Let us discuss now the gauged NJL model in curved spacetime. 

The assumption that not elementary ones but composite bound states are mainly relevant in quantum 
cosmology motivated the study of the chiral symmetry breaking under the influence of the external 
gravitational field (see previous sections.) Such an investigation has been done via explicit calculation 
of the effective potential for the composite field ('tpip). Our purpose here will be to discuss the same 
question on RG language for the gauged NJL model using its equivalence with the gauged Higgs- Yukawa 
model (5-1) (the explicit loop calculations in such a model are extremely hard to do). 

We start from the gauged NJL model with four-fermion coupling constant G in curved spacetime 



-G 2 

4 » v 



N f n f 

jyivfD^ + g Y^^i 

i=l i=l 



(176) 



As usually one introduces an auxiliary field a to replace the NJL model by the equivalent Higgs- Yukawa 
model. As is shown in Ref. pf|, it is possible within the RG approach to use a set of boundary conditions 
for the effective couplings of the gauge Higgs- Yukawa model at t\ = ln(A/jUo) (where A is the UV-cut 
off of the gauged NJL model) in order to identify the gauged NJL model with the gauge Higgs- Yukawa 
model. Taking into account Eqs. (5.8) the explicit expressions for these compositeness conditions (which 
are specified in Ref. p5|) are given as (see also Eqs. (4-7) and (4-8) of Ref. p4j): 
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(177) 



where t < t\. Hence, the system (5T) with cut-off A is equivalent to the system (5T9) with the same 
cut-off when the running coupling constants are given by (5-20). In this sense the gauged NJL model 
may be called a renormalizable one. In addition, as it has been shown in Ref. |Q, one should have also 
a compositeness condition for the mass, whose most convenient form will be (for b ~ > +0 p4[): 



m 2 (t) 



2a 



y\{t) f 



1 



54(A) 



(178) 



where G = ((47r) 2 /a)g4(A)/A 2 , 34(A) is a dimensionless constant, and w = 1 — a/(2ac), a = ao and 
where y\(t) is given by the first of Eqs. (5-20) at b — > +0. 

The above compositeness conditions define the gauged NJL model as the gauge Higgs- Yukawa model 
in flat spacetime. Now, since we are working in curved spacetime (in linear curvature approximation), 
one has to add the compositeness condition for £(£a). Making the calculations in the same way as in 
Ref. HI] one obtains 

£(*a) = \ ■ (179) 

Obviously, the compositeness condition for £(t) is again the same one as in the non-gauged NJL model 
[52]| (for finite corrections see also Ref. J53|). Hence, the conformal invariance due to the renormalization 
group ]5(| takes place again. Analyzing the RG equation (5T3) for £(t) we will find that in all cases 
(i.e., for the general situation, where y(t) is given by (5-20), for the case with 6-^0 and with the 
corresponding y(t), and also for the case ao — > 0, i.e., the non-gauged NJL model) we have 



(180) 
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Now, one can study the possibility of chiral symmetry breaking in the model under discussion using 
the effective potential language, or more precisely, the RG-improved effective potential (see Refs. |E3 
and |5j| and references therein). Because the technique to study the RG-improved effective potential 
is widely known for the flat space |5f| as well as for curved spacetimc ji9[ |5lJ we will not give any 
details of its derivation. 

Using the fact that the effective potential satisfies the RG equation, one can explicitly solve this 
equation by the method of characteristics, and we find: 

V(g, y, A, m 2 , t a, M ) = V(g(t),y(t), A(t), m 2 (i), £(<), a(t), pe*) , (181) 

where the effective coupling constants g(t), a(t) are defined by the RG equations (5-5), (5-13) (at the 
scale fie*) and corresponding RG equations for m 2 (t), a(t) written in Ref. |Q for the case of the above 
gauge-Higgs- Yukawa model (t is left unspecified for the moment). As boundary condition it is convenient 
to use the one-loop effective potential (5-18). 

In the case of the gauged NJL model one has to substitute into the effective potential (5-18) the 
effective couplings fulfilling compositeness conditions in accordance with Eq. (5-24). In this way, one 
obtains the RG improved effective potential in the gauged NJL model. In flat spacetime such calculation 
has been already done in Ref. [Q for the case of fixed gauge coupling (b — > +0), so we may use for the 
running coupling constants (except scalar-gravitational coupling constant) the results which are known 
from flat spacetime calculations. 

Taking into account that the condition 

A*e* = Ai F (t), (182) 
serves actually for finding t, we get for the case of fixed gauge coupling 

/ / * x l/(2-io) 

eJ = fejd) , (183) 

where ^f(m) = f^F^t = 0). Using (5.20), (5-23) and (5-26) in the case 6^0 one gets (for the case of flat 
space see expression (6-13) of Ref. Q): 
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24/i 2 \ fi J \ a 

where g\ = w. Thus, we have got the RG-improved effective potential for the gauged NJL model (with 
finite cut-off) in curved spacetime. 

Taking the limit A — > oo we will get the renormalized effective potential of the gauged NJL model in 
curved spacetime: 

(4tt) 2 V _ If 1 1 \ yl a 2 ( M ) 
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(185) 



where y 2 = (47r) 2 a/2aac and the four-fermion coupling renormalization has been done; in particular 
g\ R is a finite constant (for discussion of that renormalization see Refs. [Q and ]5q|). 

Note that in flat space (i? = 0) the potential (5-28) in the same approach has been obtained in 
Ref. (5(| (for non-renormalized potential see also Ref. |5^| ) using the ladder SD equation. Hence, the 
RG-improved effective potential may serve as a very useful tool to study the non-perturbative effects on 
equal footing with the SD equation. It is really surprising that RG improved effective potential gives 
the same results as ladder SD equation. 

Using the finite effective potential (5-28) we are able to discuss the chiral symmetry breaking in the 
gauged NJL model under consideration. In flat spacetime, the possibility of chiral symmetry breaking is 
defined by the sign of the first term in (5-28). That gives the value of the critical four-fermion coupling 
constant 54(A) = w. In curved spacetime the situation is more complicated as one has additional terms 
in the effective potential (5-28), even in the linear curvature approximation. 

Taking into account that 

1 3C 2 (R) 2 1 3 N c g 2 

an 2n 4n 
a 3 N c g A 

and introducing x = er(/i)//i, one can rewrite the quadratic part of the potential (5-28) as follows: 



(4tt) 2 fl / 1 1 \ R a c 



x 1 . (187) 



Relation (5-30) determines the way to estimate the chiral symmetry breaking for the gauged NJL model 
in curved spacetime. 

In particular, in flat spacetime and for (cut-off) dependent effective potential we observe that the 
chiral symmetry is broken for 

' 1 < or — l — - -L < . (188) 



54(A) g\ 54fl(M) g\ R 

Hence, the critical value of the four-fermion coupling constant which divides the chiral symmetric and 
non-symmetric phase is g\ = w. 

In curved spacetime, chiral symmetry is always broken if 

1 1 \ an R 

' ~ -o < • (189) 



guiM g\ R ) 6a fi' 

When this condition (5-32) is valid one can easily find the curvature-induced dynamical fermion mass. 
In a similar way one can find the chiral symmetry breaking condition for cut-off dependent effective 
potential, 



fi 2 J V.94(A) g\) 6(i 2 a 



i/ac 



< . (190) 



Thus, we have got the condition for chiral symmetry breaking in terms of the dimensionless curvature 
R = Rj fi 2 . From (5-32) we see that the critical coupling constant depends on the curvature. For negative 
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curvature the critical coupling is higher and one has a greater chance to find the system in the phase 
with broken chiral symmetry. As an example of spaces being in correspondence with (5-33) one can 
consider the inflationary universe S 4 with small curvature. 

Note that the critical value of the curvature at which symmetry breaking is absent is defined according 
to (5-32) by 

i?cr 6a 

,2 ~ 



(191) 



M a cg AR 

At all positive curvatures below Rqt as well as at small negative curvatures the chiral symmetry is 
broken. 

It is quite remarkable that such a simple quantum condition of chiral symmetry breaking in gauged 
NJL model in curved spacetime is obtained explicitly. Up to now such a simple symmetry breaking 
tree-level condition in curved spacetime has been known only for the Higgs sector of (5-11): 



£R + m 2 
A 



(192) 



where £i? + m 2 < 0. 

Thus we discussed the gauged NJL model in curved spacetime using quite standard RG language. 
The effective potential for composite fermions is found explicitly and its phase structure is discussed. 
Some analytical results are obtained explicitly (in particularly, the condition of chiral symmetry breaking 
in the case of fixed gauge coupling). There are different possibilities to extend above methods. First, 
one can consider other gauged NJL models (with more scalars) as equivalent to gauge-Higgs- Yukawa 
models typical for GUTs. Second, one can consider other gravitational backgrounds (say, de Sitter 
space). There is no problem also to generalize this approach for the situation when external gravity is 
sufficiently strong (see Ref. |p8| ). 



5.2 Higher derivative four-fermion model in curved spacetime 

As it has been mentioned several times, four-fermion models are considered now as non-renormalizablc 
effective theories where the presence of an ultraviolet cut-off A at loop diagrams is a necessary condition. 
Then there are different possibilities to extend the NJL model. Among such possibilities, a quite inter- 
esting one is the introduction in the original Lagrangian of higher derivative terms in the four-fermion 
interaction []59| or the kinetic term. |6(| It maybe shown |H| that the physics of such higher derivative 
four-fermion model is still equivalent to the physics of SM. Note also that the inclusion of higher deriva- 
tives in effective theories gives the possibility to take into account the structural effects of the medium 
and external fields. In the present subsection we discuss the effective potential and phase structure of 
higher derivative four-fermion model in curved spacetime. We mainly follow Ref. [p2| . 

We start by presenting the model which we set out to study (it was introduced in Ref. |59|]). Its 
Lagrangian in curved spacetime is given by 



C = ■0i7 M (x)V Al -0 + 
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(193) 



where N c is the number of fermionic species, A a cut-off parameter, and Ai and A2 are coupling constants. 
We work in the 1/A C expansion scheme. 

By introducing some auxiliary fields xi and %2 we can give a description of this nonrenormalizable 
theory by means of the action 



5' 



(194) 
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Our purpose is to study the influence of external gravity on the dynamical breaking and restoration 



patterns of the symmetry possessed by the Lagrangian (193), which are given by the transformations 
Tp — ► 75-0 and ip — > ^75 . If we refer to the action (194), the symmetry is realized by adding the following 
transformations for the auxiliary fields: X x ~^ — Xi an d X2 — * ~X2- 
The effective potential of this model in the N c — ► 00 limit is given by 

Ai A2 

with 

\-i ■ 



y = A 2 I X2 + X2 
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(195) 



Vi = i (Vol)' 



InDct (i W ~ m+ 2V3 X 2^j 



where Vol is the volume of the spacetime and m — xi + V^X2 ■ Thus, to leading order in the 1/iV, 
expansion we have the gap equations, as follows 

dVi 
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(197) 



i f -m + 2V3x2fi 2 /A 2 ' 
First of all, we have to calculate the fermionic propagator, which satisfies 

i JJ -m + 2V3 X 2j2^ Q(x, x') = 5(x, x') . 

Once this propagator is obtained (apart from terms which disappear under the tr operation) it will be 
immediate to produce Eqs. ( |196[ ) and (197) in explicit form and also the dependence of the effective 
potential itself on the fields xi an d X2 ■ The details of this rather lengthy derivation arc written in the 
Appendix of Ref . |62| , where we give details about the calculation of the effective potential up to terms 
linear in the curvature. The technique is actually equivalent to weak curvature expansion method given 
in §3. (The difference is in the presence of higher derivatives in propagator equation). 

The final outcome is shown below in natural variables, which are obtained by using the following 
definitions 

R Xi X2 . /r V 

r = -p ) x\ = -T- , x 2 = ^-, a = x x + V3a;2 , V = . 

The 'dimensionless' effective potential v maybe calculated as ]62| 



(198) 



h 



a 

r 2 



2 

a r 



32tt 2 



96tt 2 2V6 h 



V3^ 

47T 2 



2vr 2 



(199) 



where l\ — A1/A 2 



h 2vr 2 
h = A 2 /A 2 . 
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Note that in deriving this expression we have taken into account only 
terms up to quartic order on the fields a and X2 ■ The gap equations — which one obtains by differentiating 
v with respect to x\ and x 2 and equating the results to zero — are, respectively 
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One can study now the influence of gravity on the symmetry breaking pattern of the theory. A simple 
inspection reveals that a positive curvature tends to protect the symmetry of the vacuum, and also a 
negative one may trigger the breaking of the symmetry. 

In flat spacetime the symmetry is broken whenever either Ai or A2 are greater than 87r 2 /A 2 , that is 
why we shall take in the following k\ and to be l\ — 871" 2 and I2 — 8ir 2 respectively. We may illustrate 
these remarks by studying the evolution of the minimum, given by two coordinates (which we choose to 
be a and X2) in several circumstances. Before including gravity, it is worth discussing in more detail the 
situation in flat space. In fact, several cases may be analyzed (see Ref. [[39|). To introduce this study it 
is convenient to define 

8n 2 " 



A' 1 



where A^ means A,: > 8n 2 (for the + sign) or Ai < 8n 2 (for the — sign). One sticks here to the case 
when fii <§; A. We only repeat the two situations given when fj, 2 > 3fi 2 and (A^ , A^), or /1 2 < 3fi 2 and 
(A^, A^). In both cases the results may be summarized by 
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(202) 



ln(A 2 /r 



As a first example, consider the case in which the coupling constants are such that the symmetry is 
not broken in flat space, that is to say, A^ < 8tt 2 . One can see here the situation is modified as we move 
from negative to positive curvature by the numerical analysis. ]62| We observe that there is a negative 
value of the curvature above which the symmetry is restored (phase transition.) 

Note that usual four-fermion model was studied under the influence of a gravitational field in different 
situations (see §3). The conclusion of this section was always that if the coupling constant is greater 
than the critical value in flat space, there is a first order phase transition at some positive value of the 
curvature. In the present case there is no actual contradiction. The point is that, till now, we have 
concentrated ourselves in cases where the coupling constants are around 8tt 2 . We observe that if we 
explore regions in the space of parameters where A2 is much smaller, and Ai is greater than 8ir 2 — which 
would correspond to a limit where our theory approaches the Gross-Neveu model — then there is a first 
order phase transition for some positive value of the curvature. Similarly, one can study numerically 
phase structure of the theory for other values of parameters and curvature. 

There is an interesting question about this model — namely whether it would be possible for some 
generalization of it to be represented in renormalizable form. For A2 = it is possible to transform 
theory to Yukawa-type model which describes near the critical point the physics of chiral symmetry 



breaking. Note that the higher derivative term in (193) acts against such generalization at A2 =/= 



Finally, note that in the same way, one can study other higher derivative four-fermion models in 
curved spacetime. 

5.3 Supersymmetric NJL model in curved spacetime 

Next we consider the supersymmetric NJL model non-minimally interacting with external supergravity. 
It is known that in flat spacetime |33| the fermion and the boson loop effects in this model are canceled 
out. Thus the supersymmetric NJL model does not show the dynamical chiral symmetry breaking in 
flat spacetime. Hence, could we expect that external (super) gravitational field could be the source of 
phase transitions in SUSY NJL model? 
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Below, we will discuss the curvature induced crural symmetry breaking in supersymmetric NJL model 
following closely Ref. |64[ ]. 

We will start with the action of SUSY NJL model in an external supergravitational background. This 
model can be considered as local supersymmetry generalization of the model given in Ref. |63| : 



S = / d*zE- 1 [QQ + Q c Q c + X (Q c Q)(QQ c ) + i 1 {QQ) + UQ c Q c )] , 



(203) 



Q c a carry the color index a = 1, • • • , N and belong to the representations of 



M , u M is supertetrade, £1 and £2 are non- minimal coupling constants of SUSY NJL 



where chiral superfields Q 
SU(N), E = BevE*, 

model with external supergravity. We follow the notation of book |65| . Note also that there exists more 
general form of non-minimal interaction above SUSY theory with supergravity but we consider only a 
simplest variant. 

After the standard introduction of auxiliary superfields in a similar way which is explained in §2 and 



rewriting the action (203) in component fields we will limit ourselves to purely gravitational background. 
In the leading order of the 1/iV-expansion the spinor auxiliary fields are dropped away since they may 
contribute only to next-to- leading order terms in the l/7V-expansion. Then the action to start with takes 
the form: 



S = 
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+V(«7 M V A 
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(VV p 



(204) 



where p 2 = a 2 + ir 2 is an auxiliary scalar as in the original NJL model, ip is N component Dirac spinor, 
£1 = (1 + £i)/6, £2 = (1 + ^2)/6. The minimal interaction with external supergravity corresponds to 
£1 = £2 = 1/6. The start point action has the chiral symmetry. If the auxiliary field p develops the 
non-vanishing vacuum expectation value, (p) = m ^ 0, the fcrmion ip and the scalar <j) acquire the 
dynamical mass m and the chiral symmetry is eventually bro ken] 1 " 4 ] 

To find the phase structure of the model given by action (204) we introduce an effective potential. 
To evaluate the effective potential ones start with generating functional of Green functions. 
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An internal line of the p-propagator has no contribution in the leading order of the 1/iV-expansion. 
Assuming that the p is slowly variating field and applying the 1/iV-expansion method the effective 
potential for p is found to be 



V(p) = 



1 2 

— p -Itt 
AO 



ds S(x, x; s) 



-2i 



sds [Gi(x,x; s) + G 2 (x 1 x; s)] + O [ — 



(206) 



14 We started with the locally supersymmetric action (203) in order to justify the consideration of the Lagrangia n (B04 ) 
and clarify its origin. However, after restriction to pure gravitational background the theory with Lagrangian (E04I) is 
not supersymmetric since the gravitino field ( fermionic superpartner for gravitational field ) is absent. However, taking 
into account that we started from supersymmetric theory (in flat spacetime) we continue to call our model as SUSY NJL 
model. 
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where S(x, X] s) and Gi(x,x] s) are the spinor and scalar two-point functions respectively. It should be 
noted that the effective potential (206) is normalized as V(0) = 0. 

We will evaluate the effective potential taking into account the terms up to linear curvature and using 
local momentum representation of propagators (see §3 above). Hence the effective potential reads 
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Thus for £i + $2 = 1/2 the fermion loop contribution is cancelled with the boson loop contribution. To 
evaluate the integration in Eq. (207) we perform the Wick rotation p° — > ip° 
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sds 
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Applying the Schwinger proper time method [fl2| / is rewritten as 
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Since the integration over t is divergent around t ~ 0, we introduce the proper time cut-off A and find 

R 

- I 

1/A 2 

p 2 Ei 



I 



2(4tt) 2 
R 



p_ 

A 2 



A 2 (e~ 



2(4tt) 2 

where Ei(— x) is the exponential-integral function which is defined by 
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Hence the effective potential in the leading order of the 1 /TV-expansion for the supersymmetric NJL 
model in curved spacetime reads 



where is 



p 2 Ei(-^) + A 2 (e^-1) 
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(212) 
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The ground state of the system corresponds to the minimum of the effective potential. To find the 
ground state we evaluate the effective potential with varying the f(£,i,£,2)R^o- |@ 
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Figure 19: Behaviour of the dynamically generated mass m with varying the &)RXq. 



In flat spacetime R — the effective potential ( [21 2| ) takes a very simple form of quadratic function. 
It is evident that the minimum of the effective potential will be only at p = for arbitrary values of the 
nonminimal coupling parameters. 



For /(£i,£2)i?Ao > taking into account the property (211) of the exponential-integral function we 
see that the effective potenial (212) is non- negative and takes minimum at p = 0. 

To understand the phase structure of the SUSY NJL model for a negative /(£i, £,2)RXq we will analyse 
the effective potential more precisely. The dynamically generated mass m of fermion tjj and scalar (f> is 
given by t he v alue of p at the minimum of the effective potential. Stationary condition for the effective 
potential (212) is given by 
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Thus the dynamical mass m = (p) satisfies 

16vr 2 
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(214) 
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This equation has a solution only for a negative &)RXo- In Fig- [jJI the dynamically generated 
mass is plotted as a function of £,2)RXq. The dynamical mass m — (p) which corresponds to the 
order parameter smoothly disappears as the /(£i> £2)i?Ao increases. 

The critical curvature i?cr which divides the symmetric phase and the broken phase is given by the 
curvature where the order parameter m disappears. If we take the limit rn/A — > in Eq. (215), we find 
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(216) 



Due to the stability of the ground state the coupling Ao must be real and positive. Hence the critical 
curvature is given by Rex = and the chiral symmetry is broken down for an arbitrary negative /(£i, £,2)R 
and an arbitrary positive Aq . 
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For m/A <C 1 the dynamical mass m is obtained by 



m 2 



= A 



167r \ „ / m 



cxp( /(a,6)i?A -^ +o U 



,2 



(217) 



The dynamical mass is suppressed exponentially for a small negative ^-R^o- 
Thus it is found that the chiral symmetry is broken down for an arbitrary negative /(£i, £2)-/? and an 
arbitrary positive Ao in SUSY NJL model within our approximation. For £1 +£2 < 1/2 the broken phase 
is realised in a spacetime with an arbitrary negative curvature and the symmetric phase is realized for 
an arbitrary positive curvature. On the other hand there is only the broken phase in positive curvature 
spacetime and the symmetric phase at a negative curvature for £1 + £2 > 1/2. In both cases the critical 
curvature is given by i?cr = 0. The dynamically generated fermion mass exponentially disappears at the 
critical point. If the dynamically generated fermion mass is extremely small, the large R 2 /m 2 compared 
with R may spoil the validity of the linear curvature approximation. The terms of the order 0(R 2 /m 2 ) 
may change the exponential behavior of the dynamical mass near the critical point. 

Note finally that there are few possibilities to generalize the results of the present consideration. First 
of all, it could be interesting to study phase structure of SUSY NJL model on non-trivial supergravita- 
tional background (with non-zero gravitino). Second, it could be interesting to investigate the gauged 
SUSY NJL model where compositeness conditions a la Bardeen-Hill-Lindner jl5| maybe implemented. 
That could lead to the formulation of compositeness condition for £1 and £2 which presumably should 
be given by asymptotically (super) conformal invariant values. |t], |5^] 



6 Thermal and curvature induced phase transition 

At the early stage of the universe the combined effects of the temperature and curvature to the DSB 
may play an important role. In the present sections we will apply the similar analysis as in §2 at finite 
temperature and curvature. To show the combined effects of the temperature and curvature to the 
DSB, we fix the coupling constant A above the critical one and see whether the broken chiral symmetry 
is restored in an environment of the high temperature and/or large curvature]^] To introduce the 
temperature in the theory it is supposed that the system is in equilibrium. This assumption is not 
accepted in a general curved spacetime. In the spacetime which has no time evolution the equilibrium 
state can be defined. We then restrict ourselves in the positive curvature spacetime R ® S ^ 1 and the 
negative curvature spacetime R^H ^ 1 . To find the ground state at finite temperature and curvature we 
calculate the effective potential and analyze its stationary condition by the gap equation in the leading 
order of the 1/iV expansion. 

15 In Ref. |66| the combined effects of the temperature and curvature to the DSB is discussed in the positive weak 
curvature limit. In the paper the variable k is introduced by using the trace formula: 



(4tts) d / 2 J (2tt) 



Inserting Eq. ( 21£ ) into Eq. (|25|) the effective potential is rewritten as 



1 r 1 f d D k 2 

V{a) = — a 2 - itr / ds / — -e~ sk S(x,x;s). (219) 

2A J Q (Aits)"/ 2 J {2tv) u 

The temperature T and chemical potential /i are introduced by the replacement 

fc° - y/{uj n - i^) 2 , I |£ -» k B Tj2 , (220) 



where kg is the Boltzmann constant and u> n = (2n + l)ivkgT. The variable k in Eq. (219) does not always correspond 
to the momentum of the fermion field in a general curved spacetime. Thus above method is useful only in flat spacetime 
or special spacetime where the variable k is regarded as the momentum of the fermion field. So we do not use the trace 
formula to introduce temperature here. 
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In this section we mainly follow Ref. 
dynamical symmetry breaking. 



and discuss the thermal and curvature effects to the 



6.1 Effective potential at finite temperature 

It is known that the broken chiral symmetry is restored for a sufficiently high temperature through 
the second order phase transition in flat spacetime.Jlf], Here we introduce the effect of the finite 
temperature to the four-fermion models and give the explicit expression of the effective potential at 
finite temperature. 

As we have seen in §2, the effective potential is expressed by the two-point function S(x, x; s) of a 
massive free fermion. The two-point function at finite temperature is defined by 
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where E a is the energy in the state specified by quantum number a respectively, (3 — 1/UbT with ks 
the Boltzmann constant and T the temperature. 

Following the standard procedure of the dealing with Matsubara Green function, the two-point func- 
tion at finite temperature is obtained from the one at T — by the Wick rotation and the replacements 
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In Minkowski space the effective potential at finite temperature in the leading order of the 1/N expansion 
reads M 
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If we perform the integration over k 7 we obtain 
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Performing a summation and integrating over angle variables and s in Eq. (|223|), we get 
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Comparing Eq. (224) with Eq. (|225| ) we find the following relation. 
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This relation will be used for numerical calculation of the effective potential at finite temperature in 
curved spacetime. 

Starting from the theory with the broken chiral symmetry at vanishing T and evaluating the effective 
potential at finite temperature, it is known that the broken chiral symmetry is restored at a critical 
temperature Xcr through the second order phase transition. The critical temperature is given by p0[ 
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Below, we investigate the combined effects of the temperature and curvature on R®S and R®H 
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6.2 Positive curvature space (R <g> S ) 

In the positive curvature space, R ® S ^ 1 , the effective potential is given by (4-60). According to the 
definition of the two-point function at finite temperature (221), we obtain the effective potential in the 
space at finite temperature by the replacements (222) 
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where a n is defined by 
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Evaluating the effective potential (228) we will find the phase structure of the model at finite temperature 
in positive curvature space. 

For numerical calculations we need the finit e expressi on o f the effective potential in summation and 
integration. Inserting Eq. (2-33) and Eq. ( 226 ) into Eq.(228) the renormalized effective potential reads 
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(230) 



In this representation of the effective potential the divergence is canceled out in the summation. 

The phase structure of the theory is obtained by observing the minimum of the effective potential. 
The necessary condition for the minimum of the effective potential is given by the gap equation: 
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If the gap equation has a non-trivial solution which corresponds to the minimum of the effective potential, 
the chiral symmetry is broken down and the dynamical fermion mass is generated. The non-trivial 
solution of the gap equation is given by 
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Acr is defined in Eq. (^) and m corresponds to the dynamically generated 



where a n = a^J m 2 + 
fermion mass. 

In Ref. |67| the typical behaviors of the dynamical fermion mass m was evaluated numerically as a 
function of temperature T or curvature R. It was found that the broken chiral symmetry is restored for 
a sufficiently high temperature and large curvature and only the second order phase transition is realized 
with varying temperature and/or curvature for 2 < D < 4. 

Since the dynamical fermion mass smoothly disappears at the critical point for second order phase 
transition, the critical line on T-R plane is given by the massless limit of Eq. (232). To find the equation 
for the critical line in an analytic form we take the limit m — » in Eq. (232) and find 
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where C(z,a) is the generalized zeta function, u>cr n = (2ji + l)7r//3cr, acr n = a cr l w cr,J and mo is the 
dynamical fermion mass in Minkowski spacetime at T = 0. The critical lines are shown in Fig. 

The two-dimensional spacetime R <g> S 1 is a flat compact spacetime, R — 0. Thus the symmetry 
restoration which is caused by decreasing a is induced by the finite size effect of the compact space. In 
four dimensions the effective potential at T = is divergent. However the thermal effect gives only a 
finite correction to the effective potential. Thus the critical temperature Tcr at R = is divergent at the 
four-dimensional limit. Near four dimensions the curvature effect seems to give the main contribution 
to the phase transition. But it results from the non-renormalizability of the four-fermion model. In a 
renormalizable theory the situation must be changed. 



6.3 Negative curvature space (R (g) H D l ) 

In a negative curvature spacetime the chiral symmetry is always broken down irrespective of A at T = 0. 
Can the thermal effect restore the symmetry in the negative curvature spacetime? Here we consider 
the model at finite temperature in the negative curvature spacetime, R® H ^ 1 , for both A > Acr and 
A < Acr- The manifold R ® H 0-1 is defined by the metric 

ds 2 = dt 2 + a 2 (d6 2 + sinh 2 6 dCl D _ 2 ) . (233) 
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Figure 20: The phase diagram at D = 2.0, 2.5, 3.0, 3.5. 



It is a constant curvature spacetime with negative curvature 

R = -(D -!)(£>- 2)a~ 2 . 



(234) 



respectively (2 < D < 4). 

According to the similar method used in the previous subsection the effective potential at finite 



temperature in R (g> H D 1 reads 



trl 



3-D °° 



o (47r)(^- 1 )/ 2 /3 



E 



■ D-l 



'3-D 



3-D 



(235) 



R <X> i? 1 is equivalent to the two-dimensional Minkowski space R 2 . At the two-dimensional limit of 
Eq. (235) the effective potential in two-dimensional Minkowski space is reproduced. Because of the 
convenience for numerical calculations we rewrite the effective potential ( |235|) in the same form described 
in the previous subsection. Inserting Eq. (2-33) and Eq. ([226|) into Eq. (|235|) we get 
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To find the minimum of the effective potential (236) we analyze the non-trivial solution of the gap 
equation. Substituting Eq. (236) to Eq. (231) the gap equation reads 
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Figure 21: The phase diagram at D=2.5, 3.0, 3.5 for A > A c 
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Evaluating the gap equation (237) the dynamical fermion mass m is obtained. The dynamical 
fermion mass smoothly disappears as the temperature increases with the curvature fixed for both A > Acr 
and A < Acr- Then the broken chiral symmetry is restored for a sufficiently high temperature. On 
the other hand the dynamical fermion mass becomes heavier as the curvature \R\ decreases with the 
temperature fixed. Calculating Eq. ( |227| ) in three dimensions the critical temperature for a — > oo is 
given by 

*bT ci . = ^- 2 . (238) 

The curvature effects enhance the symmetry breaking on H 2 . Hence there is only the broken phase 
for k B T < 1/(2 ln2) in the model A > A CT - For k B T > 1/(2 ln2) or A < A C r the dynamical fermion 
mass is smoothly generated as the curvature \R\ increases and the chiral symmetry is broken down by 
the curvature effect. After the same analysis done in arbitrary dimensions 2 < D < 4 we find the same 
behavior for the dynamical fermion mass. The thermal effect restores the broken chiral symmetry while 
the negative curvature effect breaks the chiral symmetry. Only the second order phase transition occurs 
with varying temperature and curvature in 2 < D < 4. In two dimensions Eq. (237) has the same 
behavior in Minkowski space. 

For the second order phase transition the critical poin t is obtained by the massless limit of the gap 
equation. Taking the massless limit m — > in Eq. (237) we find the equation that gives the relation 
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Figure 22: The phase diagram at D = 2.5,3.0,3.5 for A < Acr 



between critical temperature /?cr = l/(fcsTcr) and critical radius acr- 
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(239) 



acr i 



«cr„r + a C r„) 



= 0. 



Evaluating Eq. (239) numerically we draw the phase diagram of the four-fermion model with varying 
temperature and/or curvature on fl® H 13-1 at D = 2.5, 3.0, 3.5 in Figs, ^l] and In those figures the 
normalization scale m is taken to the value defined in Eq. (Et]) for A > Acr and 



mo = /i 



(4tt 



.D/2 



1 



,l/(D-2) 



trir(l - D/2) X R 



-D + l 



(240) 



for A < Acr- At a scale T ~ a -1 ~ too the thermal effect gives the main contribution to the phase 
structure for D < 3 and the symmetric phase is realized. At four-dimensional limit Tcr is divergent for 
A > Acr- It comes from the divergence which appears at the four-dimensional limit. Thus it may be a 
result caused mainly by the non-renormalizability of the theory. 

Thus we have investigated the phase structure of the four-fermion model at finite temperature and 
curvature in arbitrary dimensions (2 < D < 4). In positive curvature space the curvature effect restores 
the broken chiral symmetry. On the contrary the curvature effect enhances the chiral symmetry breaking 
in negative curvature space. At finite temperature the lower limit of the momentum for the fermion field 
(k° < it/ (3) appears from the antiperiodicity and then the long range effect is suppressed. Thus the 
thermal effect restore the broken symmetry even in a negative curvature spacetime. 



Gl 



In two spacetime dimensions it is impossible to introduce the combined effect of the temperature 
and curvature within our method because of the assumption of the equilibrium. This assumption is not 
fulfilled in inflationary expanding universe, but we may discuss the phase transition at early universe 
using our results. 



7 Finite size and topological effects 

In the previous sections we discussed four-fermion models in curved spacetime. The main motivation 
to do such a study was the cosmological applications. However, one may expect that the very early 
universe had a non-trivial topology (the introduction of the temperature in the hot universe maybe 
classified as the similar type of effects). This renders interesting to investigate the phenomenon of 
dynamical chiral symmetry breaking in curved spacetime with non-trivial topology. This section will be 
devoted to such a study on the background M ' 1 ® S 1 where S 1 is the one-dimensional sphere and 
.M D_1 is a (D — l)-dimensional arbitrary curved manifold of trivial topology (see Ref. [p9| ). 



7.1 



S 1 



Effective potential in M. D 

We start again from the action of four-fermion models in curved spacetime defined in (Q) or (^) . Making 
the explicit calculations in the same way as in §2 we easily get the effective potential in the leading order 
of 1/iV expansion: 

Using the known expression for the propagator of a free Dirac field in a weakly-varying gravitational 



background (see expression (67)), we will then write the effective potential with accuracy up to linear 
curvature terms. In the spacetime M. ^ 1 ® S 1 there is no constraint for the boundary condition |7(]. [nj 
along the compactified direction. Some of the independent spin structures are allowed in our universe. 
If we take the antiperiodic boundary condition along the compactified direction, the field theory in 
M 13 ^ 1 ® S 1 is equivalent to the finite temperature field theory (see §6). Here we consider the fermion 
fields with periodic and antiperiodic boundary conditions independently. Thus, the effective potential is 
found to be 



Via) = — a 2 -itr / sds- V 
2A T. 
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7U-12 
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(241) 



where one should integrate over kg, • • • , ku-2 and sum over the coordinate feo-i, which is given by: 
&d_i = (2n + S P: i) tt/L. 



Evaluating the effective potential (241) we study the phase structure of the four-fermion models in 
weakly curved spacetime with nontrivial topology below. 



7.2 Flat spacetime with nontrivial topology (R D 1 <g> S 1 ) 

Now it is not so clear if the symmetry restoration is mainly caused by the curvature effect or finite 
size effect in an arbitrary dimension. In flat spacetime with nontrivial topology we clearly see only the 
finite size effect. In Refs. [72, [73|, [74| the dynamical symmetry breaking has been investigated in flat 
spacetime with nontrivial topology. In the present subsection we calculate the effective potential ( ^4l| ) 
in flat spacetime R 0-1 ® S 1 to see only a finite size effect. 



G2 



In R D 1 (g> S 1 the effective potential (241) reduces to 

If we perform the integration in Eq. ( |242| ) we find that (TgJ 

- 2 trl 



(242) 



2A 2(4tt)( Z5 - 1 )/ 2 
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(243) 



Thus the effective potential is modified b y fin ite size effect which is interpreted as a generalized version 
of the Casimir effect. |75| Analyzing Eq. ( |243| ) it is found that the divergences in the effective potential 
are of the same form in Minkowski spacetime. Thus we need not renormalize the parts which include the 
finite size effects contributions. We apply the same renormalization condition as that shown in Eq. ( |32"| ) 
and obtain the renormalized effective potential as in Minkowski spacetime. The renormalized effective 



potential is given by replacing Ao in Eq. (243) with the renormalized one A in Eq. (£ 
If we apply the expression (|243[) in the gap equation (pCl), we find that 
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In Eq. (244) we neglected the trivial solution m = 0. 

In the case of antiperiodic boundary condition it is well-known that the field theory in i? D_1 



(244) 



>S l is 



equivalent to the finite temperature field theory. The effective potential ( 244 ) is in agreement with that 
of the finite temperature four-fermion theory (see for example Ref. |l6|| and references there in) with 
recourse to the relation between the size L of the compactified direction and the temperature T 
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(245) 



with kg the Boltzmann constant. Evaluating the effective potential (243) and the gap equation (E4 



one can study numerically phase structure of the theory and find that the second order phase transition 
takes place and the broken chiral symmetry is restored for sufficiently small L < Lqt at the large N limit. 



It agrees with the known result at finite temperature. Analyzing the gap equation (244) the critical size 
is given by 

'2r((3-£>)/2) lVP-2) 



L cr mo = 2tt 



{2 6 ~ u - 1)C(3 - D) 



(246) 



ySFT(l - D/2) ' 

where mo is the dynamical fermion mass in Minkowski spacetime. It is equal to the well-known formula 
of the critical temperature for the four-fermion model. (lq] In Fig. 23 the critical size is plotted as a 



function of dimension D for A > Acr by the use of Eq. (246). 

On the other hand finite size L has an opposite effect for the fermion field with periodic boundary 



condition. In this case the first derivative of the effective potential (244) has a negative value for D = 2 
at the limit a — > 

dV trl 1 „ . „„, 

— ->-- __< 0. (247) 
da a ^ Q 2 L 

For 2 < D < 3 the first derivative of V vanishes but the second derivative of V has a negative value at 
the limit a — > 
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(248) 
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Figure 23: Critical size L C r as the function of D for A > Acr in the case of the anti-periodic boundary 
condition. 
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Figure 24: Critical size L CT as the function of D for A < A C r in the case of the periodic boundary 
condition. 



64 



Thus the symmetric phase cannot be realized for 2 < D < 3 even if the coupling constant A is sufficiently 
small. 

For 3 < D < 4 the first derivative of V also vanishes at the limit a — > but the sign of the second 



derivative of V depends on the coupling constant A. Evaluating the effective potential ( 243 ) numerically 
it is observed that the chiral symmetry is broken down for a sufficiently small length L in the case 
A < Acr, while only the broken phase is realized for A > Acr- 



Analyzing the gap equation (244) the critical value of L is given by 

l/(D-2) 



Lcim' = 2ir 



2r((3-D)/2) . 
V^T(l-D/2) a3 °\ 



(249) 



where m' is defined in Eq. (p4|). In Fig. 24 we present the critical size Lqt as a function of dimension 
D by the use of Eq. ( |249| ) . In drawing Fig. 24 we fix the coupling constant A < A cr . Since the finite 
size effect gives a finite correction to the effective potential even at the four-dimensional limit, the finite 
size effects disappear at D — » 4 as is shown in Figs. 23 and 24. 

Thus the boundary condition of the fermion field has large contributions to the phase structure. If 
the fermion field possesses the antiperiodic boundary condition the finite size effect raise the vacuum 
energy V((<r)). Hence the broken chiral symmetry is restored for a sufficiently small L. The fermion field 
which possesses the periodic boundary condition has an opposite finite size effect. The vacuum energy 
V((a)) reduces by the finite size effect. Thus the chiral symmetry is broken down for a small L. 

We can interpret above results in the following way. Under the antiperiodic boundary condition 
minimum momentum allowed for the fermion field becomes p = n/L. Then the infrared cut-off induced 
in a finite L spacctimc. Since the lower momentum modes have an essential role to break the chiral 
symmetry, the symmetry restoration occurs for a sufficiently small L. On the other hand the vanishing 
momentum is allowed under the periodic boundary condition. In £g) S 1 a fermion field ip{x) is able 

to interact with ?p(x + nL) where n is an arbitrary integer. Summing up all the correlations, {ip(x)ip(x)) , 
{^}{x)tjj{x + L)) , (<ip(x)'ip(x + 2L)), ■ ■ •, the vacuum expectation value of the composite operator increases 
as L decreases. It is understood as a dimensional reduction. Compactifying one direction to the size L 
in D-dimensional space, it looks (D— I)-dimensional space for particles with Compton wavelength much 
larger than the size L. In the lower dimensional space the influence from the lower momentum fermion 
exceeds. Then the finite size effect breaks the chiral symmetry under the periodic boundary condition. 

7.3 Curved spacetime with nontrivial topology (Ai 3 (g> S 1 ) 

Now we are interested in the combined effect of the finite size and curvature. The fermion field with an 
antiperiodic boundary condition has been already discussed on compact spaces S D and R(g> S " 1 in §4. 
Thus we consider the universe where both the independent spin structures exist below. 



For this purpose we evaluate the effective potential (241) in weakly curved spacetime. Since we 
are interested in the theory in four dimensions, we restrict ourselves in four spacetime dimensions and 
calculate the effective potential by using the cut-off regularization. Notice also that a summation over 
the two inequivalent spin structures which are admitted by the spacetime A4 3 x S 1 will be performed. 
Of course, one can consider also the corrections corresponding to periodic and antiperiodic (non-zero 
temperature) boundary conditions independently. 



Integration over s in Eq. (241) is immediate. To perform the momentum integration, one first makes 
the Wick rotation (A; = ifc 4 ) and puts then a cut-off to regularize the resulting expressions. In our case, 
we simply restrict 

(fc 4 ) 2 + (fc 1 ) 2 + (k 2 ) 2 < A 2 , 

so that our cut-off is different, when compared with the cut-off for the case of trivial topology discussed 
in §3.2. 

From now on, we shall call V± the contribution to the effective potential which comes from the 



logarithm of the determinant of the operator which appears in Eq. (241). After carrying out the 
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integrations over s and the momenta, we are led to the following expression for the contribution to 
V\ coming from the p — case, namely, purely periodic boundary conditions (this corresponds to the 
contribution to V(a) obtained by taking only the p = term in Eq. ( |24l| )). Of course, we want to study 
V\ , which is given by the sum over the two values of p, but — as we discuss below — V\ may be written 
down immediately once the p = contribution has been worked out. We obtain 
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(250) 



To simplify this expression we use standard techniques drawn from complex analysis, such as the 
expression 
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One has just to identify the function / for each term in Eq. (|250j) and then perform the integrals in 
(251). Another remark is in order here: when computing the term 
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it is better to rewrite it as 
lim 
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since now the expression within square brackets satisfies the properties which justify the use of Eq. 

As we said before, once V[ has been computed, one may write immediately Vf 1 and Vi, which is 
given (by definition) by Vi = V[~° + Vf~ , because 
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and, then 
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Now it is apparent that the physics displayed by the model in the case of purely periodic boundary 
conditions and in the one where we consider both spin structures will be essentially the same, since 



vr L ° + v?= L i = 2v^i 



Thus we see that both cases are related by a trivial rescaling of the length and an overall factor which 
multiples V\. 

However, the last remark does not apply to the case of purely antiperiodic boundary conditions, where 
we only take p = 1 — which gives the thermodynamics of a system in three-dimensional space with trivial 
topology. In fact, from the last expression we get 
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Henceforth, we shall concentrate below only on the analysis of the case in which both spin structures 
are taken into account, as they appear in Eq. (241). 

The first term on the right-hand side of Eq. ( |25l| ) may be computed without difficulty in all cases. 
The second term is, in general, rather more involved. In order to simplify this contribution as much as 
possible in the cases which come from Eq. ( |250| ), one should pay careful attention to the determination 
of the integrand along the contour of integration. After some work, the final result is found to be 
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where x = a/A, I = LA and r = R/A 2 . 

The value of the field a which satisfies the gap equation, V'(a) — dV(a)/da = 0, gives a dynamical 
mass to the fermions. This last equation, when written in terms of the natural variables x, I, r and c 
(c = AA 2 ), reads 







A 4 



x 
~2c 



1 



3vr 2 + x 2 1 - exp (-2ZVT 



arcsinh [ — J — \J\ + 



1 : 
2^ 7f 



67 



4 

3^ 



X 



1 



3 (2<T 

2rx 
3(27r) 2 



VI + x 2 exp (2JV1 + Z 2 ) - 1 
arcsinh'^ 1 

1 



- 3x 



/' T 2_ 2 ,2\ 1 /2 



exp(2Zr) - 1 



xj VI + x 2 
1 



4(3tt) 2 (1 + x: 2 ) 3/2 



Vl + x 2 exp (2/Vl + x 2 ) - 1 
dr 1 



2(3tt)' 
21 



Vr 2 - x 2 exp(2Zr) - 1 
a; 1 



(255) 



(1 + ^2)3/2 exp (2/VTT^) -1 
exp (2/Vl + x 2 ) 



1 + x 2 



(exp (2ZyT 



Let us discuss some limiting case of above expression (]255|) . Assuming LA <C 1, let us expand V/A 4 
in powers of I. Assuming now that + x 2 -C 1, we readily obtain 
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where g = XqA /I. Using (256) one can study the phase structure of the model. 



7.3.1 Phase structure in Ad c 



S 1 



Because we are making expansion on the curvature, i.e., r -C 1, it is worthwhile studying first the case 
r = 0. Here it can be proven that it is impossible to have a first order phase transition, in fact there is 
a second order phase transition at gcr — vr 2 /2. For g > gcr there is chiral symmetry breaking and the 
symmetry is restored for g < gcr\ in other words, for a given value of A, the symmetry is restored when I 
grows beyond a critical point. In the broken phase the order parameter is given by x^ r = n (2/tt 2 — 1/g) 
or, in terms of the generated mass mgen = 2A/n — Ltt/Xq. It is straightforward to study the behavior 
of the value of the effective potential at x Dr , for varying g or I, in this limit of small compactification 
length. One finds 
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Now we shall study the influence that the presence of curvature has on this behavior. It is immediate 
to notice that for negative values of the curvature, the symmetry is always broken (the slope of the 
effective potential is r/24-7T at the origin.) Moreover, for fixed positive values of r (which is kept always 
small), as I grows there is a first-order phase transition: it has actually lost its continuous character. 
Now the critical value of the parameter g is 

7T 2 / 27T [7 

9cT = Y\ Tv 8 
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The approximation is consistent with the small-curvature limit. The difference between the values of 
the order parameter x in the broken and disordered phases at the phase transition is given by x = r/8 
(retaining again only the first correction coming from the curvature). 

We can now, on the other hand, study the situation when g is held fixed and the curvature takes on 
different values. 

If g < ir 2 /2 there is a continuous phase transition at r — (the symmetry is broken for negative 
values of r and is restored for positive curvature). As r approaches from below, the order parameter 
tends to zero according to the expression 

I r I 



24ttA ' 

or rngen = A _1 |i?|/247rA, being A = l/g — 2/ir 2 . In deriving this result we assume that | r |<C 1 and 
that | r |< A. 

One may also consider the case A < or, equivalently, g > ir 2 /2. To be consistent with the require- 
ment of small curvature near the critical point, we have to assume now that | A |<C 1. In this setting 
one may expect that, as r grows, there will be a phase transition from the ordered to the disordered 
phase at some positive value of the curvature. Retaining again the first terms of the expansions only, we 
obtain 



(3tt) 2 a2 
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We shall here consider again the phase structure in the fixed-curvature case but corresponding now 
to the opposite limit, when L is large. 



It is easy to see from Eq. (255) that, dropping exponentially vanishing contributions, one may 
approximate it with 
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where c is, as before, c = AoA 2 . 

From the analysis of gap equation it is trivial to see that the symmetry is broken when 

Otherwise, the symmetry is respected by the vacuum. 

If one takes into account the presence of a background gravitational field, one can check that the 
effective potential at large L is becoming 
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To analyze phase structure let us consider A < 7r 2 /Ag. Then the order parameter maybe found approxi- 
mately as 

^(jTj^/c-i)-^ 

(remember that this is valid for negative r). The symmetry is restored for positive values of r. The 
numerical analysis for this and some other cases maybe found in Ref. p9[. 
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Thus we discussed phase structure of four-fermion model on M. <£> S . In compact flat spacctimc 
R ^ 1 ® S 1 it was found that the boundary condition for the fermion field drastically changes the phase 
structure. In the case of an antiperiodic boundary condition finite size effect restores the broken chiral 
symmetry, while the finite size effect breaks the chiral symmetry for the fermion field with a periodic 
boundary condition. In A4 3 (£)S 1 we found the possibility of (topology combined with curvature)-induced 
phase transitions. 



8 Chiral symmetry breaking in curved spacetime with magnetic 
field 

The external electromagnetic (EM) fields play an important role as a possible probe of quantum field 
theory (for a review of quantum field theory in an external EM fields, see Ref. [Q). It is quite well- 
known fact that an external magnetic field supports chiral symmetry breaking in four-fermion models, 
(see Refs. |77| ~ and references there in). There was some activity on the study of phase structure 
of four-fermion models in an external magnetic field (for a review, see for example Ref. |80| ). 

From another point it maybe possible that strong primordial magnetic fields should be considered 
on equal footing with the strong curvature in the early universe. Then, different consequences of the 
combined effect of the magnetic and gravitational fields may occur (for a example, it maybe significant 
increase in the number of created particles in the early universe filled with the constant EM field 0]). 
Hence, for cosmological applications it could be interesting to discuss the chiral symmetry breaking 
phenomenon in four-fermion models under the action of external gravitational and magnetic fields. This 
section will be devoted to such situation where gravitational field will be taken in linear curvature 
approach (see §3). We follow Refs. Q and Q. 



8.1 Three-dimensional four-fermion model with magnetic field 

We start from three-dimensional four-fermion model. The effective potential maybe easily written as 
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where at first step there is no magnetic field yet. It is convenient to work in terms of the derivative from 
V(<j). Then, using the local momentum representation of the propagator (see §3), we get 
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Making use of Wick rotation and calculating the trace in (|260|), we obtain 
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In derivation of Eq. ( 261 ) we have used the expression for proper-time representation 
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and after that, the ultraviolet proper-time cut-off A 2 has been introduced in the low limit of proper-time 
integral. This cut-off is different from the cut-off used in §3. 

Performing the integration over s and over a, we get (compare with Ref. p3[ or §3 where other 
regularization has been used) : 
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where erfc(x) = (2/y / 7r) 



dt. Hence, we obtained the effective potential in three-dimensional NJL 



model in curved spacetime using proper-time cut-off. 

The renormalizcd effective potential may be found in the limit A 



as following:]^ 
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Let us analyze now the phase structur e of the potential (264). In the absence of curvature, R = 0, 
for A > the minimum of the potential ( 264 ) is given by a = 0. Hence, there is no chiral symmetry 
breaking. For A < the chiral symmetry is broken, the dynamically generated fermion mass is given as 



m 



7T 

A 



(266) 



When the curvature is not zero we find the following picture. Let first A > 0. Then the ground state 
(and dynamically generated mass) is defined by 



7T 

2A 



7T Z 

A 2 



R 
~6 



(267) 



One can see that unlike the case of flat space for positive A we have the chiral symmetry breaking. We 
also see the possibility of curvature-induced phase transitions. The critical curvature is given by i?cr = 
(flat space). For negative curvature, R < i?cr = 0, we observe the chiral symmetry breaking, while for 
positive curvatures symmetry is not broken. 

For negative four-fermion coupling constant A < we get the following ground state: 



2A 



A 2 



R 
~6 



(268) 



The critical curvature is defined by the condition: Rci — 97r 2 /2A 2 = (9/2)m§. Between < R < i?cr 
the chiral symmetry is broken. 

Now, let us discuss the situation when four-fermion model is considered in curved spacetime with 
external magnetic field. That means that spinor covariant derivative also contains the electromagnetic 
piece. Treating the magnetic field exactly, (8^] one can easily find the effective potential for our model. 

16 Three-dimensional four-fermion model is renormalizable in the sense of 1/jV expansion. Jl^] Thus the theory is defined 
independent of the regularization method. If we perform the finite renormalization 

1 1 2 

A ~* A ~ 7T ' 

and put trl = 4, the effective potential (jrfj) is reproduced. 
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Considering now four-fermion model in curved spacctimc with magnetic field, we again work in linear 
curvature approximation as above (but making no approximations for an external magnetic field). More- 
over, we take into account only leading contribution on curvature which does not depend on magnetic 
field, i.e., the contribution discussed above. (One can show that for not a very large magnetic field the 
curvature correction depending explicitly from magnetic field is not essential). Then, using Eq. ( |264j ) 
and the results of the calculation of three-dimensional four-fermion effective potential in an external 
magnetic field, |77], one can get: 



2A 3tt 24tt 4tt 3 /2 /„ s 3/2 
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(269) 



where H is magnetic field and e is elect ric ch arge. 

The renormalized effective potential ( ^69|) may be represented as following 



2A 24tt 2vr 2ir ^ \ 2' 2eH 

where properties of generalized zeta- function may be found in Ref. |$0| . 

Working in frames of our approximation and considering also eH — ► 0, we find that for positive A 

A f eH R\ , , 

^~A-2-2i) ' (271) 

Hence, for R — chiral symmetry breaking due to magnetic field occurs. Negative curvature increases 
chiral symmetry breaking. However, positive curvature acts against the chiral symmetry breaking. On 
the critical line of phase diagram 

— ~eH (272) 
12 v ; 

the restoration of chiral symmetry breaking occurs. Hence, in this example magnetic field and gravita- 
tional field act in the opposite directions with respect to chiral symmetry breaking. Gravity tends to 
restore the chiral symmetry while magnetic field tends to break it. On the critical line both effects are 
compensated and there is no chiral symmetry breaking. 

Similarly, one can analyze other choices for A, eH (see Ref. [g2| for more details). 

8.2 Four-dimensional NJL model in curved spacetime with magnetic field 

Let us start again from the action for the NJL model in curved spacetime: 

S = J d^x^—g Ui-y»(x)V^ + ^ [{Uf + [&1^Y] } ■ (8 ■ 15) 

Standard calculation gives the effective potential as (where as usually it is enough to put pseudo-scalar 
7r = as is shown in §2.) 

2 

V(a) = ^—+i \nV>et[i^(x)V^ -a]. (8 • 16) 

2Ao 

Then calculation similar to the one in previous subsection gives the effective potential in proper-time 
cut-off 

" 1 R 
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where Ei(— x) is defined in Eq.(5-54). 

Expanding Eq. (8-17) and keeping only terms which are not zero at A — > oo we get 
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2A (4tt) 2 
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Thus, we have got the effective potential with proper-time cut-off. 
Using Eq. (8-18) the gap equation is found as follows 
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This gap equation defines the possibility of chiral symmetry breaking in curved spacetime (in linear 
curvature approximation). If our weakly curved spacetime is filled by the magnetic field (the covariant 



derivative is now V, 



-Bx2§fj,i) one can get the following effective potential 
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-HlnDet 
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In the absence of the gravitational field (R = 0) , the effective potential corresponds to flat space situation 
where the magnetic field is treated exactly. [[12] In the absence of the magnetic field (B = 0) we are back 
to the potential (8-17). In addition, in the linear curvature terms the effect of the magnetic field is not 
taken into account as in the previous subsection. 

Making the calculation of the integrals in Eq. (8-20) up to 0(1/A 2 ), and taking the derivative with 
respect to a one gets the gap equation, dV/du = 0, as follows 
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Using this gap equation one can study the dynamical symmetry breaking in different cases. In some 
cases it can be given analytically, for example, for values of the coupling constant A much below the 
critical value, i.e., 



A < 



4^ 
A 2 ' 



One can find (supposing that the second term on the r.h.s. of (8-21) is the leading one) 
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and finally, the dynamically generated fermionic mass is given by 



R , A 2 
In — 



u 



22) 



23) 



\eB\~ 

4:7T 



-\eB\/A 2 



1 

A 2 



-_R/(12A 2 ) 



exp 



4tt 2 
X^A 2 



l/[|eB|/A 2 +_R/(12A 2 )] 



24) 



73 



In the absence of the magnetic field it gives the analytic expression for the dynamically generated 
fermionic mass due to the curvature 
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One can see that positive curvature tends to make the first term in (8-24) less, i.e., it acts against the 
dynamical symmetry breaking. At the same time, the negative curvature always favors the dynamical 
chiral symmetry breaking in accordance with the explicit calculations in as external gravitational field 
(see §3). 

The numerical analysis shows the similar qualitative behavior (for more details, see Ref. |p2[ ). Note 
finally that in the same way one consider more complicated situations when both gravitational and 
magnetic fields are treated exactly (see, for example, Ref. [[^5] where two dimensional Gross-Neveu 
model on sphere with magnetic monopole has been considered). 



9 Quantum gravity and four-fermion models 

In the present section we give few examples where four-fermion models maybe relevant to study the 
dynamics of quantum gravity. In particular we apply the effective potential of D = 2 four-fermion model 
to investigate D — 2 black hole solution of dilatonic gravity. We discuss conformal factor dynamics in 
frame of 1/N expansion. The effective potential for composite gravitino on de Sitter background is also 
investigated. 



9.1 Semiclassical approach in the dilatonic gravity with Gross-Neveu model 

Motivated by the idea that 2D quantum gravity is easier to study than AD quantum gravity there was 
recentely much of interest in 2D quantum gravity and its black holes solutions. |8(| |7|, [8£| |i| |9(J We 
consider now one of such solutions in the 2D dilatonic gravity with Gross-Neveu model following Ref. 
@. 

The action to start with is 
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-gd x 
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fc=i 
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(273) 



where </> is the dilaton field and Gross-Neveu action is the same as Eq. (2-12) 
Notice that if we redefine the metric by the relation 

_ -2<j> - 

g^iv — e g^iu, 



the gravitational part of the action (272) takes the form 



(274) 



St 
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-gd 2 xe- 2 ' t ' (i? + Ag^d^dA + Ae" 2 *) , 



(275) 



which is similar to the CGHS action |g8| (the cosmological co nsta nt term is a fixed number in Ref. |p8f 
while it is accompanied by the Liouville type potential in Eq. (275).) However, we should emphasize that 
the theory with action (|273| ) is different from the theory with action (|275| ) coupled to the Gross-Neveu 
fermions (with metric tensor g^u) already at the classical level. The difference becomes more serious in 
the semiclassical approach. 

Let us start our discussion of the theory described by Eq. (273) in the semiclassical approach (l/N 
expansion where dilatonic quantum effects are subleading) . We will assume that background spinors are 
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absent. In the semiclassical approach we have to add to the action (|273j) the conformal anomaly term 

MM 



s c = 



i 
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—gd x 
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(276) 



where Z is the auxiliary scalar field. For the semiclassical argument we need the G ross- Neveu effective 
potential V(a, R) in curved spacetime which repla ces classical potential a 2 /2 in Eq. (272). The dilatonic 
equation of motion derived from the action ( |273| ) 

i? + A = 0, (277) 



apparently constrains the curvature to be constant. Using the Gross-Neveu effective potential in curved 
spacetime with constant curvature in the leading order of the 1/N expansion, our effective action in the 
semiclassical approach reads 



S cS=^J V=9~d 2 x [e~ 2 *(R + A) - V(a, R)] + S c 
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where V{a, R) is given by Eq. (3-15) with cr-independent part to be discared. 
Applying the conformal gauge 
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(279) 



we derive the following semiclassical equations of motion from the effective action (278) (with constant 
curvature according to Eq. (277) and without spinors): 
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(280) 



N 



2d 2 



+ ^(d-Z) 2 + 



{d± - 2d- P d-)Z = 0, 



(281) 



and 
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dV(a,-A) 
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The general solution of Eq. (283) is well known 
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where F± is an arbitrary holomorphic (anti-holomorphic) function. The conformal gauge in Eq. ( ^86| ) 
can be fixed completely by choosing = x ± . 
The solution of Eq. (|284|) is given by 



P + U + (x + ) + U-(x ), 



(287) 



where u± is an arbitrary holomorphic (anti-holomorphic) function. Using Eqs. (283) and (284) in Eq. 
( p82| ), we get 

d+d-e- 2 ^ + 4e 2p (e~ 2 ^ - a) = 0, (288) 



where 
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(289) 



and Co is a solution of Eq. (285). Note that the quantity Aa plays the role of the effective cosmological 
constant. It may be interesting to note that in the weak curvature limit the effective cosmological 
constant disappears (a = 0) if 
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(290) 



Accordingly we observe that the quantum effects in the Gross-Neveu model provides us with the way to 
solve the cosmological constant problem within the model under investigation. 
The general solutiuon of Eq. (28S) is 
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Here the function u± and v± are determined by solving the equations for T ++ and T The equation 

for T ++ with the choice F± = x ± reads 



x+ 
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(292) 



and the same for T with the change of index + to 

We rewrite these two equations in the form 
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(293) 



where we omit indices + and — , p(x) is an unknown function and { , } is the Schwarzian derivative. By 
solving Eq. (|293]) for v and u and substituting the solutions into Eq. (291), we finally have 



-2rj>(x + ,x ) 



N 
9G 



C+x+ + 



C 



i A + - 
1 — — x^x 



+- 



8" : 

'■V 



dy_ 

y 



{p+{z),z}dz- 



a 



z 2 {p + (z), z}dz 



1 f v 

-- J z 2 {p + (z),z}dz 



ctv + — + y { P +(z),z}dz 



2+(+--), 



(294) 



76 



where Cf and Cf are arbitrary constants. For the simplest case in which p±(x ± ) = x ± , the solution 
for dilaton field is given by 
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(295) 



This is the main result of this subsection. We found the static solution for our model in an explicit form. 
In order to compare it with the CGHS-Witten black hole solution, (8^] we make the redefinition (274) of 
the metric. The new spacetime interval squared becomes 
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This new solution may be regarded as a modified version of the CGHS-Witten black hole, 
the simplest case where all arbitrary constants are chosen to vanish [eT 2 ^ = a), we get 



ds 2 



-dx + dx 



A , 
a 1 + —xx 



(296) 

(297) 

(298) 
| For 

(299) 



It is interesting to note that for A = in Eq. (273) the semiclassical solution becomes 

p = 0, 

2 2 f 27r 
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(300) 



where Vb = V(ao, R = 0) and Co is an arbitrary constant. In terms of the new metric (274) we obtain 

dx + dx~ 



ds 2 



Co - ^V x + x~ 



(301) 



Equation ([301]) exactly coincides with the black hole solution of Refs. [|86| and [ p8| . 

Thus we presented the example of 2D black hole solution with constant curvature in dilatonic gravity 
with Gross-Neveu model. It is interesting to remark that one can use the renormalization group improved 
effective potential instead of one-loop Gross-Neveu effective potential. Then, one has only change the 
coupling constant (A = Ng 2 ) in the semiclassical effective action by the running coupling constant 
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where t is the renormalization group parameter. Working in the regime of strong curvature, t = 
In R/ p 2 . |92[] Hence, asymptotic freedom in Gross-Neveu model is induced by black hole curvature. 
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Stronger black hole curvature induces stronger asymptotic freedom as in d = 4 gauge theories in curved 
spacetime.0 In other words, in vicinity of constant curvature black hole, Gross-Neveu model tends to 
become free theory. 

Note also that one can consider other types of dilatonic gravity interacting with Gross-Neveu model 
and investigate the constant curvature black holes solutions there. However, to study other types of 
black hole solutions one should calculate the effective potential of Gross-Neveu model in an arbitrary 
curved spacetime. 



9.2 Conformal factor dynamics in l/N expansion 



The conformal factor dynamics describes the conformally-flat solution of the gravitational theories at 
classical level. At the quantum level the dynamics of conformal factor (induced by the conformal 
anomaly) was suggested as the tool for the description of quantum gravity in the infrared phase, pi 
Conformal factor dynamics gives rise to the effective potential for the conformal factor which seems to be 
quite relevant in quantum gravity. j94j In particular it maybe responsible for infrared phase transitions. It 
is very interesting to note that composite bound states which are typical for four-fermion theories maybe 
also the origin of the conformal factor dynamics. We follow Ref. |)5| in discussion of this question. 
Our starting point is the two-dimensional theory with action 
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d 2 x- 
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i\> {i-i^{x)V^ -m)ip + R 
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where the massive A-component spinor ip is considered to be a quantum field. The gravitational field, 
on the other hand, may be either classical or quantum. We also consider the conformal parametrization 
of the metric 



(304) 



where p is the conformal factor (in general it is p = p{x)) and 77 M „ is the flat fiducial metric. In the QG 
case, the choice (304) corresponds to the gauge fixing. Substituting (304) into (303) one gets, at the 
classical level, 

S = J d 2 x x (Wm - mp ) x -±p 2 , (305) 
where \ — p 1 / 2 ^. Rescaling p — > A 1 / 2 /}, we get 



S = 



(306) 



where h = mk 1//2 . As one can see, action (|306| ) has the form that is typical for the Gross-Neveu model 
(p ~ xXi see Eq. (2-12)). The dynamics of this model are quite well known: asymptotic freedom in the 
UV limit 



where t is the RG parameter. However, since (306) describes also the dynamics of the conformal factor, 
the interpretation of the function h 2 (t) is now completely different from the original interpretation. The 
h 2 (t) here is a combination of the fermionic mass in (303) and of the two-dimensional cosmological 
constant A. Using the anomalous scaling dimension one gets the running composite field 



p(x,t) = p(x) (1 + h 2 Nt/n) 



-1/2 



(308) 



The conformal factor has acquired the classical dimension after the rescaling p — * A 1 / 2 p. Hence, one may 
argue that the t dependence is due completely to that of the cosmological (dimensional) constant, i.e., 
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A(i) ~ A (1 + t^Nt/Tr)^ 1 . Therefore, there appears to be a screening of the cosmological constant in the 
1 /N expansion in the UV regime (quantum gravitational corrections are subleading in 1 /N expantion) . 

One can also investigate specific features of the effective potential for the conformal factor, which 
coincides with the GN effective potential.^) In particular, the appearance of a minimum, i.e., a non-zero 
vacuum expectation value (v.e.v.), for the conformal factor 
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is interesting. After having shown the possibility to study the dynamics of the conformal factor in two 
dimensions as the dynamics of the GN model, we now turn to the four-dimensional theory which is 
physically more interesting. 

We start from the multiplicatively renormalizable theory j?j with action 



S = 



xv-g 



A 



R 



w 
a7 



UR 2 
3Ai 



(310) 



where ip is an A^-component spinor and W the square of the Weyl tensor. The gravitational field may 
be chosen to be classical or quantum, and the theory remains multiplicatively renormalizable in both 
cases. 



We shall work again with the conformal parametrization (304) for the four-dimensional metric. In the 
case of four-dimensional QG this does not fix the gauge, contrary to what happens in two dimensions, 
but it can still be considered as a convenient background. Q Rewriting action (|310D, we get 
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where x = p 3 ^ 2 4> an d cr = lnp. In this way we have got the classical theory for the conformal factor. 



At the quantum level, the theory (311) may be considered as a n effective theory for QG (see also Ref. 



|93J). If we drop the p-terms with derivatives from action (311), we obtain a model that is reminiscent 



of the NJL model (where, of course, owing to the absence of the M 2 p 2 -term, it is p ~ (xx) 1 ^ 3 ) 



Now we are going to study the theory (311) in the large- A limit, while concentrating our attention on 
the RG and low-derivative terms in (311). The higher-derivative terms are actually of lesser importance, 
moreover, they simply disappear in the subsequent analysis of the effective potential for the conformal 
factor. First of all, we rescale p — > \f\2 p/n and denote h = m/e/vT2 and A = Ak 2 /6. Thus, 
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By integrating over the fermionic field, we get the effective potential for the conformal factor at large A^: 



V(p) = iNTrln (i^d» - hp) , 



(313) 



where p is constant. Supposing that, as usually, for large A^ A scales as A ~ AA^, where A does not 
depend on N, and using a finite cut-off p (see also §3) we get (notice that N has been factored out) 
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The v.e.v. of the conformal factor can be found from (314) as the solution of the equation 
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One can present a sample of numerical values of the solution p 2 //i 2 (which lie between and 1 for 
0.04644 < A < 1). Notice that from the point of view of the original theory, a non-zero v.e.v. for p is 



more acceptable physically, because for p = the conformal parametrization (304) becomes degenerate. 



Now we turn to the study of the renormalization structure of theory (312), which is rather non-trivial. 
There are a few different ways to renormalize this theory, the actual problem being the fact that p is 
dimensionless. 



We may consider a theory (312) — as it stands — which will correspond eventually to some QG phase. 



After renormalization of p (taking into account the negative sign for {dp) 2 ), we obtain 

4Nh 4 _ 3A 2 - &N\h 2 - 48Nh 
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The theory has now a UV stable fixed point (i — > +oo), where the behavior of h 2 and A is 



& = -773J. & = • ( 316 ) 



h 2 (t) ~ 'A X(t) ~ (317) 
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Hence, we now obtain a decrease of the cosmological constant in the UV limit. Notice that in comparing 
this theory with the standard scalar self-interacting theory, we do not have here physical restrictions on 
the sign of A, and a negative sign is perfectly acceptable. That indicates to the possibility of solution of 
cosmological constant problem in frames of conformal factor dynamics based on four-fermion like theory. 

9.3 Effective action in iV=l supergravity 

In §5.3 we already discussed supersymmetric extention of NJL model. There are also other possibilities 
to introduce supersymmetry into the theory. In particular, in present subsection we will discuss N = 1 
supergravity on de Sitter background. We calculate the effective action for non-zero gravitino condensate 
o ~ ( - 0aL ab- 0t>) ■ Hence, such a model maybe considered also as kind of four-fermion (four-gravitino 
model). Note, however, that unlike the above discussion we work in one-loop approximation (not in 
1 /TV-expansion) . We follow Ref. |96| below. 

The Lagrangian for the theory of N = 1 supergravity is 

C = ™iZ(e, ^) + \^ vpa i>^l^ P i>o + \{Al -S 2 - P 2 ) , (318) 

where T^ v = ify'*, y u ], i?(e, ip) = i?(e) + iifc 4 (?/' /j r MI/ -;/v) 2 + - ■ • (in this expression the terms corresponding 
to the interaction of the gravitino with the gravitational field and four-fermion terms with insertions 
containing the 75 matrix are not written in the explicit form), Va = Vx(e), A^, S, P is the minimal set 



of auxiliary fields. Let us represent the Lagrangian (318) in the form:|96l p7\ 



L = ~^2 R ( e > ^ + ^" CT VV757,V p ^ + a 2 

-VUk*@ ll T>»'1> v ) + J(A 2 - S 2 - P 2 ) , (319) 



wher e a is the auxiliary scalar field. Let us calculate the effective action in the theory with the Lagrangian 

it scalar field a. If the effective equations 

= 0, (320) 



( 319 ) in de Sitter spacetime (R = 4A) with constant scalar field a. If the effective equations 

ST 
5a 



§ = 0, (321) 

have the real solutions a 7^ 0, then these solutions indicate dynamical supersymmetry breaking. |97| 
Equation ( |321| ) represents the effective equation for the gravitational field which depends on the param- 
eter A only. 
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Let us write g^ v — > g^ + h^ v ,a — > a + a q , where g^ is the de Sitter space metric, a =const, and 
h^ u ,a q are quantum fields. The background values of other fields of the theory we put equal to zero. 
One can obtain that the sector of quantum fields (h^, a q ) does not interact with the gravitino sector. 

At first one can integrate over the field a q in the functional integral. Then the additional term 
A£ CTg = — <7 2 /i 2 /4 appears in the expansion of the classical gravitational Lagrangian on the quantum 
field hfi U . 

In order to calculate the effective action, one needs to add a gauge-fixing term to the Lagrangian. We 
shall choose the one-parameter gauge: 
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where 7 is the gauge parameter, 7=1 corresponds to the Dc Donder gauge, 7 — * 00 corresponds to the 
Landau-DeWitt gauge. Let us write the bilinear part of the Einstein Lagrangian on the S 4 background 
l9Hl with account of AC„ and C n 
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(322). 
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Here 



h = /i 



A = ifcV, 

9^K V = 0, 
VS^ = 0, 
= 0, 



(323) 

(324) 
(325) 

(326) 

(327) 
(328) 
(329) 
(330) 



Ai(X) = li(-D + A), i = 0,1, 2, 

where 1 is the unit in the space of the corresponding fields. 

Integrating over quantum fields hj^ v ,£ hi,h, taking account of the corresponding Jacobians and 
ghost contributions, we have in the limit 7 — > 00 (we omit the details of this straightforward calculation) 
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(331) 



The explicit expression for ^Sp In Aj(X), i = 0, 1, 2 is obtained in Ref. |)8| in frames of £-regularization 
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-Spin A,(A) = ~B\ In — - -(2i + 1)^(6,), 



1 



6 



(332) 
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where = 9/4 + i — 3X/A, B\ is the well-known De Witt B4 coefficient for the corresponding operator, 
F((bi) is given by Eq. (2-36) in Ref. p8| . Using (332) in (331), one can easily obtain the explicit form for 
I\ . The final expression is very complicated, so we do not write it here. It is enough for our purposes 
to have IV 



in the form (331). 

Let us now calculate the gravitino contribution to the effective action. We choose the gauge: 



(333) 



where x — V^^V- Then the Faddeev-Popov determinant is M = ID = VV + A. Represent p8f 
% - M + 7^/4, =0,^ = 4^ + (V„ - 7m V/4)C, = 0, V = 7 "V M in the bilinear part of 

the gravitino Lagrangian (31E) taking account of C^. Integrating over the fields (j>^,C,tp taking account 
of the corresponding Jacobians and ghost contribution we get in the limit 71 — > 00: 
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Sp In A 3 / 2 (to 2 



1, 



-Spin 
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; V + 2to A 1/2 



-A 



-A 1/2 (0)(V + 2m) 



SplnA 1/2 (-A), 



(334) 



where 



= — vllfccr. One can write the explicit expression for A 1 / 2 (^), A 3 / 2 (AT) using relations such 
as (332) (see Ref. Q). We do not write this complicated expression here. 



A natural question now is about the presence of the negative modes in the effective action T = 
r^„ -(- r^^. The spectra of all operators in (331) and (334) are given in Ref. p8| ]. Using the result of 
Ref. p3] for ]?/, ( [33 1[ ) one can show that the operator A2(8A/3 — k 2 a 2 ) does not include the negative 
modes if 1 < k 2 a 2 /2A < 8/3. The three last terms in ( [331| ) contain the negative modes for arbitrary 
values a and A (A > 0). Therefore, T^ contains an imaginary part for all a and A. Then T also 
contains an imaginary part. (Remember that the structure of the fermion sector does not influence the 
bosonic sector.) Thus the effective equations also contain an imaginary part. It indicates a vacuum 
instability of the considered background and the absence of dynamical supersymmetry breaking at the 
1-loop level. 

It is interesting that one can easily obtain the effective action in N = 1 supergravity in flat spacetime 
using the effective action (331) and (334). This effective action also contains an imaginary part for all 
a. 

Thus, we found the effective action for composite gravitino in N = 1 supergravity and discussed its 
properties. 

In summary, in the present section the use of four-fermion type models in the theories of quantum 
gravity has been discussed. Few particular examples presented here indicate that four-fermion models 
could find wider applications than it is usually believed. 



10 Summary and outlook 

In the present paper we discussed four-fermion models in an external gravitational field. As is well- 
known these models are extremely useful for analytical study of dynamical symmetry breaking. Our 
main purpose was to investigate the influence of the external gravitational field to dynamical symmetry 
breaking. 

Using l/A^-expansion and working in the leading order of this expansion we calculated the effective 
potential in four-fermion models in D-dimensional spacetime (2 < D < 4) for a variety of backgrounds. 
In particular, we considered the class of spacetimes where derivative expansion maybe applied. Here the 
effective potential for the composite operator ipip has been calculated taking into account terms linear 
on curvature. We also evaluated the effective potential on the background with constant curvature; 
Minkowski space, de Sitter background, anti-de Sitter background and Einstein universe. In all these 
cases an external curvature has been taken into account exactly. 
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Phase structure of the effective potential has been carefully studied using analytical and numerical 
methods. The possibility of chiral symmetry breaking due to curvature effects and curvature-induced 
phase transitions has been discussed in all detail. The phase diagrams for dynamically generated fermion 
mass have been constructed. For example, for negative curvature chiral symmetry is always broken down. 

Similar questions have been addressed for the extensions of four-fermion models; higher derivative 
four-fermion model and gauged NJL model. In particular, for gauged NJL model (which may play the 
role of SM or GUT without elementary scalars) the technique of renormalization group improvement 
has been applied. The effective potential equivalent to Schwingcr-Dyson approach has been found 
analytically. Chiral symmetry breaking has been estimated. 

We also discuss the combined influence of two effects (non-trivial topology and external curvature, 
non-vanishing magnetic field and external curvature) to phase structure of four-fermion model. The 
effective potential has again been calculated for few specific backgrounds. Phase structure has been 
analyzed. In particular we found that in some case two external effects (for example, magnetic field and 
gravitational field) may compensate each other. That leads to the restoration of chiral symmetry, which 
was originally broken. 

We presented the number of examples where an analytical study of composite states is possible even 
in the presence of non-trivial gravitational background. Having in mind the cosmological applications, 
it is straitforward now to use the effective potential obtained above for the construction of inflationary 
universe where the role of inflaton will be played by a composite state of fermion fields. Q 

In our calculation we limited ourselves to case of weakly curved or constant curvature spacetimes. In 
this case the effective action is reduced to the effective potential except for the volume factor. However, 
for more general backgrounds one has to calculate the effective action in four-fermion models. It is 
expected that instabilities (due to creation of instanton-antiinstanton condensation) may appear in such 
situation. That may completely change the phase structure of the theory. Such an extremely hard task 
as the calculation of effective action on general curved background is quite natural extension of the 
presented results. 

Another very interesting line of research is related with the study of realistic model of elementary 
particles in the early universe. One can start from some GUT scenarios which maybe reformulated 
as gauged four-fermion model in the same way as the model of §5.1. Then, the role of gravitational 
background to phase transitions maybe discussed in all detail, even in non-perturbative regime. 

Note also that methods developed in the present paper maybe applied to other theories. For example, 
in supergravities one can study the effective potential of composite gravitino (see §9.3), for gauge theories 
one can study the effective potential for composite vectors [?J and so on. 

It is interesting to note also that above approach maybe well extended to supersymmetric theories. For 
example, in §5.3 we discussed supersymmetric NJL model in curved spacetime and found the possibility 
of dynamical chiral symmetry breaking due to gravitational effects. 

It would be also extremely interesting to analyze the next-to-leading corrections to effective poten- 
tial of four-fermion model on curved background. Such multi-loop calculations may clarify the phase 
structure in many cases when the creation of multi-critical points is possible. 

Finally, it is interesting to note that in realistic situation (like early universe) the combination of 
few external effects maybe quite natural. Hence, it could be interesting to study phase structure of 
four-fermion model in curved spacetime with finite temperature , [|67| magnetic field, chemical potential 
and non-trivial topology. 
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